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Abstract—In an earlier paper [8] it was shown by visual observation of colour-change reactions 
that circulation currents within a droplet were caused by the frictional drag of the outer fluid, 
but only above a certain transition Reynolds number. The presence of circulation currents inside 
a droplet greatly accelerates mass transfer; thus, for the aqueous extraction of acetic acid from 
nitrobenzene droplets, a 2-5 to 4-5 fold increase resulted over that expected from stagnant drops 
(9). 

Observation of colour-change reaction and of flow patterns revealed by aluminium particles 
has now shown that, for droplets moving in water, a vortex becomes established in the wake of 
the drop above a Reynolds number of about 20, while the Reynolds number, at which circulation 
begins within a drop, decreases with increasing viscosity of the outer phase. Evidence is given 
supporting the earlier conclusion that droplets circulate less readily the higher the interfacial 
tension of the system. 

Trace quantities of impurities which are surface active in character may inhibit internal 
circulation, but this effect can be partially overcome by the presence of a diffusing solute readily 
miscible in both phases. 

As compared with rigid spheres of the same size and density, internal circulation may sub- 
stantially increase droplet fall velocities (an increase of about 50% was observed for one system) 
while oscillation of a prolate — oblate type will reduce the fall velocity — a reduction of at least 
12% being obtained in one case. 


Résumé—Dans un article précédent, les auteurs ont montré par l’observation visuelle les 
changements de teintes que les courants de circulation a l’intérieur d'une gouttelette sont das 
a la résistance de friction du liquide extérieur, mais seulement au-dessus d'un certain nombre 
de Reynolds de transition. La présence de ces courants internes accélére le transfert de masse : 
Par exemple, pour l'extraction de l'acide acétique par des gouttelettes de nitrobenzéne, l’augment- 
ation est de 2,5 a 4,5 fois celle que l'on attendrait de gouttes rigides. L’observation des change- 
ments de teinte et des trajectoires & l'aide de particules d’aluminium a montré qu’au-dessus 
d'un nombre de Reynolds d’environ 20, et pour des gouttelettes tombant dans l'eau, il s’établit 
un vortex dans le sillage de la goutte, tandis que le nombre de Reynolds pour lequel la circulation 
commence a l'intérieur d'une goutte décroit avec augmentation de viscosité de la phase 
extéricure. 

Il se confirme que les gouttes circulent d’autant moins facilement que la tension interfaciale 
du systéme est élevée. 


Les traces d'impuretés, de caractére tensio-actif, annulent la circulation interne, mais cet 
effet peut étre partiellement surmonté par la diffusion d'un soluté rapidement miscible dans les 
deux phases. Par comparaison avec des sphéres rigides de mémes dimensions et densités, la 
circulation interne peut accroitre de fagon substantielle la vitesse de la goutte (environ 50% dans 
un des cas) tandis qu'un mouvement oscillatoire d’allongement et d’aplatissement de la goutte 
réduira la vitesse de chute (environ 12% dans un cas). 
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THe TRANSITION FROM STREAMLINE TO 
EppyiIneg FLow arounpb Li@quIpD 
DropLets IN WATER 
When a sphere moves very slowly through a fluid, 
the flow pattern within that fluid is completely 
streamline, and the pattern in the wake is the 
same as that in front of the sphere. With increas- 
ing speed, however, a transition velocity is 
reached above which an eddy or vortex ring 
appears behind the sphere. JOHNSTONE 
Wiis [13] consider that the vortex becomes 
“unstable” at Reynolds numbers around 150, 
but this has not been confirmed by R. W. Grarron 
who worked in this field in this department. Nist 
and Porter [15] experimented with solid spheres 
in air and found that eddies first formed behind 
the sphere at a Reynolds number of 8-15. Using 
solid spheres in water, MOLLER [14] reported 
eddies forming at a Reynolds number of 300, 
though the experiments of Witiiams [18] on 
solid spheres in glycerine failed to show pro- 
duction of such eddies up to Reynolds numbers 
of 720. On the other hand, Bonn [4] states that 
the transition Reynolds number for a sphere is 

“of the order of unity.” 

We observed the transition point, using 
ethanolamine as diffusing solute passing from the 
drops into water containing about 10cc of 
standard phenolphthalein solution per litre. The 
red coloration produced by diffusion of the amine 
revealed the presence or absence of eddies behind 
the droplets. Observation was greatly facilitated 
by viewing against an illuminated background 
which was half black and half white, the droplet 
being viewed along the vertical boundary line 
between black and white, though near the tran- 
sition region the eddies were feeble and poorly 
defined. 

Two mixtures of different density and solute 
concerttrations were used (see following Table for 
composition), and for each droplet the fall velocity 
was recorded together with a comment on the 
presence or absence of eddies observed in its 
wake. Fall velocities were measured by stop 
watch and used to obtain Reynolds numbers ( Re) 
from a plot of C,/Re vs. Re, where Cp is drag 
coefficient. The method, which uses experimental 
data relating to solid spheres, was described in an 


and 


earlier paper [8], and whilst it is not valid for drops 
at higher Reynolds numbers, it is applicable here 
in view of the low Reynolds number range and 
the near-stagnancy of the drops. 





Mixture 2 


Density at 18°C 
= 1021 gm/cc 


Mizture 1 








Density at 18°C 
1-009 gm /cc 








Chloroform 144% (vol) Chloroform 25-2% (vol) 


ethanolamine 42-8°%, ethanolamine 12-0% 


amylaleohol 428% octyl alcohol 62-8% 





Mixture 1 contained a relatively high propor- 
tion of water-soluble material, ethanolamine being 
infinitely soluble and amy! alcohol 2.26% soluble 
in water. Mixture 2 contained only 12% of water 
soluble material. The drops of varying sizes were 
formed from an open tube and fell down a glass 
column 5 cm in diameter and 70 cm in length. 

The temperature of the water used was 10°C, 
its density taken as 1 gm/ce and viscosity 0-013 
poises. 

The droplets with the highest velocity in the 
region of no eddies, and with the lowest velocity 
in the eddy region, corresponded to Reynolds 
numbers (Re) on either side of the transition 
point. (Re = dup/yu, where d = drop diameter, 
v = fall velocity, p and yu are density and viscosity 
of outer phase.) 

For droplets of mixture 1, the transition 
occurred between Re of 18-6 and 22-1, and for 
droplets of mixture 2 between Re of 20-9 and 
23. Thus it appears that, for two liquids with a 
large difference in solute concentration, the tran- 
sition from streamline to eddying flow occurs 
between Reynolds numbers of 19 and 28, using 
water as continuous phase. 


Errects or Continuous Puase Viscosity 
ON TRANSITION FROM STAGNANCY TO 
CIRCULATION WITHIN DROPLETS 


Our earlier work [8] on colour change with drop- 
lets containing cobalt chloride and using water 
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40 % (vol) 





nC ,H,OH 


CHCl, 27% 

CoH, 33 %, 

CoCly 3 gm/100 ce 
+ little nCgHgOH + CHBry 
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Table 1. 
MIXTURES 


v 
(2) (3) 


As for mixture (1) 
plus extra nC,HgOH 


As for mixture (2) 
plus extra nCgH,OH 











Viscosity | Specific 
Continuous of gravity — Interfacial 
Drop phase _—_ continuous of tension 
mixture wt. %, phase continuous 16°C 
used cellofas B_ centistokes phase dynes /cm 


in water 20°C 20°C 








3-80 


' 
(1) | 0-50 11-97 10 6-24 
(1) =|) OS 24-22 10 6-24 
(1) 1-00 35-60 1-006 6-24 















Specific Viscosity 
gravity of Transition 
of disperse velocity of Transition 
disperse phase | droplet Reynolds 
phase centistokes cm /sec number Rey 


1i8-5°C 20°C 





1-045 1-61-17 64°5-71°5 








102 1-89-2-122 19-28-23 -7 





1-108 1-24 1-77-1-89 32 -— 3°59 
1-108 1-24 1-29-1-42 1-32- 1-5 
1-108 1-24 1-06-1-13 0-83- 0-91 





Nore.— Interfacial tension of mixture (1) against distilled water was 6-12 dynes/cm 


as continuous phase showed transition Reynolds 
numbers (Re,) for internal circulation ranging 
from 60 to 245. Bonp and Newron [5], on the 
other hand, in 1928, carried out experiments on 
velocities of rise and fall of droplets of mercury 
and water in media such as waterglass, syrup and 
inferred that the transition 
Spevis [17] did work on 


castor oil, and 
occurred at Re < 1. 
water in a 


oils) 


drops of glycerine and glycerine 
outer (white 
internal circulation at Re < 1. 


fairly viscous phase and 
observed Vis- 
cosities of the water-glass and syrup used by 
Bonp, and NewtTow ranged from 200-2,000 poise, 
and their castor oil had a viscosity of about 
20 poise. It was therefore decided to vary the 
viscosity of the continuous phase while attempting 
to keep other factors, such as densities of inner 
and outer phases, droplet viscosity and inter- 
facial tension, constant. 

The continuous phase comprised Cellofas B (a 
hydroxylated cellulose) dissolved in distilled 
water. The solution was centrifuged to remove a 


number of tiny hair-like fibres, and was then 
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diluted to give solutions containing 0-25, 0-5, 0-75 
and 1% Cellofas B by weight. Measurements 
showed that Cellofas B was not appreciably 
surface active and viscosities are given in Table 1. 

Droplets were formed from a mixture of 40% 
n butyl alcohol, 27% chloroform, 33% benzene 
by volume. 3gm of cobalt chloride per 100 ce 
of solution were added to give a deep blue liquid. 
More C,H,OH and a little bromoform were 
added to alter the rate of fall until smaller 
droplets were stagnant, while larger ones circu- 
lated when falling through the 1% Cellofas 
solution. This mixture was used unaltered for 
the 0.5, 0-75 and 1% Cellofas solutions, but more 
C,H,OH had to be added for the 0-25% Cellofas 
solution, to bring the drop mixture into the 
transition region. Water diffusing into the drops 
gave a blue to pink colour change, enabling intern- 
al motion to be observed. The drops were formed 
from an open tube and fell down a 60 cm column 
of about 5 cm internal diameter. 

A peculiar phenomenon was observed in which 
drops showed initial internal circulation which 
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was gradually damped out, motion ceasing first 
in the rear portion of the drop, stagnancy gradu- 
ally spreading through to the leading pole. This 
pattern was visually confirmed by the insertion 
of aluminium powder in the drops. Droplets 
absorbed no Cellofas, since their viscosities were 
unchanged after passage through the column, and 
experiments, described later, indicated that the 
explanation is some form of interfacial con- 
tamination by small quantities of surface-active 
impurity collected during passage of the drop 
down the column. 

Accordingly, observation during the first 25 em 
of fall of the occurrence of complete stagnancy 
or continual circulation was used to determine 
the transition Reynolds number from fall velo- 
cities. 

Interfacial tensions were measured by the drop 
weight method, using a thin walled glass nozzle 
of I.D. = 0-448 cm (the drop liquor wetted the 
nozzle to its inner diameter), and values were 
calculated using the correction 
method of Harkins and Brown [11)}. 
mental data appear in Table 1. 

The Reynolds number at 
stagnancy to circulation within the droplets, is 
plotted against of the continuous 
phase in Fig. 1. 

The curve appears to be a rectangular hyper- 
bola although it should be noted that the abscissa 
(viscosity) appears also in the denominator of the 
ordinate. Nevertheless, it might be expected to 
be asymtotic to the rectangular axes, since for 
zero outer phase viscosity the transition Reynolds 
number should be oo, while increasing the outer 
phase viscosity means that the frictional drag 
increases for a given velocity (in Stokesian region 
F = 8m dvp,) and hence Re, will fall. 

Confirmation of the hyperbolic nature of the 
curve is found in Fig. 2, where the log transition 
Reynolds number is plotted against the log con- 
tinuous phase viscosity. Differences in droplet 
viscosity and interfacial tension probably cause 
the values for the two lowest outer phase vis- 
cosities to lie slightly off the line, which gives the 
following approximate relationship between 
transition Reynolds number (Re,) and con- 
tinuous phase viscosity in centistokes (4) 


factors and 
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viscosity 
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EFFrects or INTERFACIAL TENSION ON 
CIRCULATION WITHIN DROPLETS 


Experiments showed that the effects of interfacial 
tension could best be detected at low fall velocities. 
Redistilled liquids with a wide range of interfacial 
tensions but with densities near that of the 
continuous phase (distilled water) were used. 

The all-glass apparatus was chemically cleaned, 
and consisted of a column 95 cm long and 4.5 «m 
in internal diameter. All drops contained grease- 
free aluminium powder to reveal internal motion. 
Drops heavier than water were formed from open 
tubes, and those lighter than water entered the 
foot of the column through various sized nozzles 
connected by cone and socket joints to a vertical 
stand pipe 130cm high and surmounted by a 
reservoir. 

The densities of indene, tetralin and 1 methy! 
naphthalene were so close to that of water that 
measurement of interfacial tensions by any of 
the normal methods proved rather difficult : 
accordingly, for qualitative indications, values 
were estimated by Antonoffs’ rule for these three 
substances. Transition velocities at which circu- 
lation began—measured by stopwatch — are 
shown in Table 2. While the values do not show 
a close correlation, they do support our carlier 
finding [8] that circulation occurs less readily the 
higher the interfacial tension of the system. 


MoDIFICATION 
TRACE 


oF CIRCULATION BY 
COMPONENTS 


A particular specimen of nitrobenzene, supplied 
as “ pure,’’ was used to form droplets containing 
aluminium particles distilled 
The 
the droplets (as 
revealed by the aluminium particles) tended to 
damp out, motion ceasing first at the rear pole 
and progressively decaying towards the leading 
pole of the drop, becoming entirely stagnant after 
115-125¢m of fall. (Droplet formation effects 
were eliminated by gently pouring the drops from 
a glass spoon held beneath the water surface.) 
This “pure” nitrobenzene was shaken with 
sodium carbonate and redistilled, whereupon 
three important changes were seen :— 


falling through 
water in a glass column 5 cm in diameter. 
circulation 


currents within 


(1) The droplets showed constant and more 
rapid internal circulation over the whole length 
of fall (143 em). 


(2) Droplets of about 0-5 cm diameter were 
about twice as oblate as those before redistilla- 
tion. This was probably a centrifugal effect due 
to greatly circulation 
the drop 


increased speed, since 
weight method 
showed the interfacial tension against water to 


be identical before and after distillation. 


determinations by 


(3) Drops of redistilled nitrobenzene about 


Table 2. 





density 
gms /cc 


viscosity 20°C 
centistokes 


interfacial tension 
dynes /cm 20°C 


transition velocity 
cms / sec 





20 
Benzyl alcohol 1043 4 


Benzaldehyde 1-046 


Indene 1-006 
Tetralin 0-973 
1 methyl naphthalene 1-025 


Cetane 0-774 





< 07 


| 


35 





* Calculated from surface tension by Antonoff's rule. 
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0-5 cm in diameter had a terminal velocity about 
25% greater than for similar size drops before 
redistillation (see later). 


The residue from the distillation was somewhat 


discoloured, but gave completely stagnant 
droplets. 

The refractive index of the redistilled 
benzene at 20-2°C was 1.5523 compared with the 
value in International Critical Tables [12] at 20°C 


of 1-5529. The refractive of the nitro- 


nitro- 


index 


benzene both before and after distillation was the 
same. 
Evidently the impurity inhibiting circulation 


is a relatively high-boiling substance, present 
only in extremely small concentration. 


INHIBITING EFFECT OF A SURFACE 
AcTIVE AGENT ON CIRCULATION 
BENAVIOUR 
A curve of drop size versus fall velocity was 
experimentally determined, and drop sizes were 
estimated from fall velocities. These drops had 


a velocity of 13-6cm/sec and a diameter of 





r 

















FALL LENGTH REQUIRED TO STOP CIRCULATION ——e 


























° 0-001 0-002 0-003 0004 0005 
wT. % AGENT IN DISTILLED WATER ——= 
Fic. 3 Disperse phase “ pure” nitrobenzene without 
redistillation, 


approximately 0-55 cm. They were formed 
beneath the surface of the water from a narrow 
glass tube, allowed to hang a moment to become 
stagnant, and then detached by a_ sideways 
jerking of the tube. This prevented any tendency 
to circulate due to formation effects. 


The disperse phase consisted of the “ pure ” 
nitrobenzene without redistillation. To the con- 
distilled water—were added 
various concentrations (0.0005% to 0.01% by 
Even 


tinuous phase 
weight) of dodecyl sodium sulphate 6. 
very small concentrations of this additive inhibited 
circulation, since the length of fall required to 
damp out circulation completely was reduced, as 
shown by Fig. 3. 


OVERCOMES EFFECTS 
AGENT 


Dirrusine So.ute 

OF A SURFACE-ACTIVE 
In these experiments the continuous phese was 
distilled water containing 0-01% by weight of 
dodecy! sodium sulphate 6. The disperse phase was 
the “‘pure”’ nitrobenzene without redistillation, and 
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FALL LENGTH REQUIRED TO STOP CIRCULATION ——= 
~ 
































C) ! 2 3 4 5 6 
VOL. % ACETIC ACID IN “PURE” NITROBENZENE —> 
Fic. 4 Continuous phase, 0-01 %, dodecyl sodium sulphate 
6 in distilled water. 
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contained various concentrations (up to 6% by 
volume) of acetic acid. The drops had a velocity 
of 12-9 em,/sec and a diameter of approximately 
0-47 cm. 

Whereas acid-free nitrobenzene droplets had 
been completely stagnant at this surface-active 
agent concentration, the presence of the acetic 
acid overcame the inhibiting effects of the agent 
and promoted circulation as shown in Fig. 4. 

Droplets of pure benzyl alcohol falling through 
distilled water, (solubility 4%) showed circulation 
(revealed by aluminium particles) almost imme- 
diately after detachment; but when they were 
formed from a glass tube previously dipped in a 
0-5% aqueous solution of sodium dodecy! sulphate 
6, circulation began only after about 30 cm of fall, 
due to the time required for dissolution of the 
surface layer of alcohol plus adsorbed film. 


STRUCTURE OF THE ADSORBED FILM 


At very low concentrations of surface-active 
material the circulating motion of the interface 
may cause the adsorbed molecules to accumulate 
at the rear of the drop surface, giving a coherent 
film which progressively extends over the inter- 
face until circulaton is stopped. 
that the rate of accumulation exceeds the rate of 


This assumes 


desorption from the film. 

The aqueous bulk concentrations of sodium 
dodecyl sulphate 6 in equilibrium with various 
concentrations at a static interface, have been 
calculated below for tentative analogy with 
conditions at the droplet interface. 

The dodecyl chain is assumed to lie flat near 
the plane of the interface and to be a cylinder [3] 
about 13-8°A in length and 20.5°A*® in cross- 


The Gibbs adsorption equation for 

C dy 
iRT dC 
is in dynes/em and the 


section [1]. 


dilute solutions is S where the 


interfacial tension, y, 
bulk concentration, C, in gm.mols/litre. Accord- 
ing to Appison [2], S approximates to the 
amount of solute in a layer one molecule deep on 


unit surface. For the concentration range 0 to 


; d 
0-5% (weight), a4 was experimentally measured 
( 


as — 6,370 in the nitrobenzene-water system, the 
Avogadro number is 6-023 x 10” and i lies 
between 1 and 2, depending upon the degree of 
ionisation of adsorbate [7] and upon the relative 
interfacial concentrations of anion and cation 
[16]. Water temperature was 20°C. 

Three types of static film structure may then 
be postulated :- 

(1) The maximum interfacial 
possible without mutual interference between 
adsorbed molecules corresponds to their forming 
the close-packed 


circles [3]. 


concentration 


diameters of hexagonally 

(2) Higher concentrations lead to increasing 
interference until a layer of tightly packed 
dodecyl chains lying flat near the interface is 
formed. 

(3) Further adsorption is accomplished by the 
two hydrophobic the molecule 
assuming an orientation more and more nearly 
normal to the interface, leading to a “* condensed 


port ions of 
se 


monolayer ”’ at saturation. 
The calculated film and bulk concentrations 
corresponding to these three types of structure 


appear in Table 3. 


Table 3. 





Film 
structure 





< 1007 x 10°10 


> 2-36 x 10° 





Interfacial concentration 


1-007 « 10°" to 2-36 x 10° > 


Equilibrium bulk concentration 


gm.mols, litre 


i=:2 





< 0-00077 
0-00077—0-00181 


> 0-00181 
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Consider the bulk concentration of 0.003% 
(weight) (0-000104 gm.mols litre). If the tran- 
sition between types 1 and 2 films prevents 
circulation, then, assuming i = 1, the bulk and 
interfacial concentrations require increasing 
0-000385 /0-000104 = 3-7 times to permit this 
film structure. If the desorption rate from the 
interface is low, the time required for adsorption 
negligible and the droplet begins circulating at 
maximum speed (i.e. renews its surface while 
falling one diameter) then, by the * accumulation 
mechanism,” this increase in interfacial con- 
centration will be achieved after travelling 3-7 
drop diameters from the nozzle (about 2 cm for a 
drop 0-55cm in diameter). The calculated 
distances of fall required to stop circulation based 
on various assumptions, are compared with the 
experimental values in Table 4 below. 

The assumptions made must lead to some 
inaccuracies, but assuming the validity of 
applying the Gibbs equation to these experiments, 
then the calculated values are probably too small 
because the drops circulate well the 
maximum possible speed [10] and because the 
time required for adsorption is appreciable in 
relation to the time of fall. 


below 


Some Factors AFFECTING FALL 
VELOCITIES IN Lieqeutp-Lieuip Systems 
Droplets were formed with a 25 cc burette fitted 
for constant head by use of a tightly fitting rubber 


bung through which passed a thin-walled glass 
tube, admitting a stream of air to a point below 
the level of the liquid in the burette. Five or six 
glass nozzles of various internal diameters were 
attached in turn to the burette with hard poly- 
thene tubing, and drop sizes were known from the 
volume of dispersed phase required to form a 
given number of drops. The drops fell through a 
glass column of 4cm internal diameter and height 
142 cm, fall velocities being measured by stop- 
watch. Three systems were studied :- 


(a) 


Redistilled * pure *’ nitrobenzene, falling in 


distilled water. 
“Pure ’”’ nitrobenzene referred to earlier, 


without redistillation, in distilled water. 


(b) 
(c) Redistilled “ pure” nitrobenzene, falling in 
0.003% (weight) solution of sodium dodecyl 
sulphate 6 in distilled water. 


The presence or absence of internal circulation 
was determined by the addition of aluminium 
particles to the drops, in parallel experiments, to 
duplicate the experimental conditions for each 
drop size. 

Discussion 

Fluid droplets differ from rigid spheres in three 
ways during free fall. 

(1) The frictional drag of the surrounding fluid 
may induce circulation of both the interface 
and interior of the droplets. 


Table 4. 





Aqueous concentration of | Calculated fall distances 
sodium dodecyl sulphate 6 | _for transition between types 


1 and 2 films at the inter- 


% (weight) 


Calculated fall distances 

for transition between types 

2 and 3 films at the inter- 
face 


Experimental distance of 
fall required to stop circu- 
lation 


cm 
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(2) The drop may depart from spherical shape, 
assuming some oblateness, due to gravity. 


(8) Prolate-oblate type of oscillation may occur. 


In Fig. 5, the experimental curves, 1, 2 and 3, 
corresponding to systems (a), (b) and (c), are 
compared with the predicted curve for rigid 
spheres of the same density, curve 4, using the 
data of Davies [6]. Shaded and unshaded points 
denote oscillating 
respectively. The experiments with aluminium 
particles showed that for curve 3 all drops were 
internally stagnant ; for curve 2 all drops showed 
weak internal circulation, while all drops of curve 


and non-oscillating drops 


1 showed vigorous internal circulation. Also, all 
except the smallest drops of curve 1 showed 
vigorous oblate-near prolate oscillation ; whereas 
all drops of curves 2 and 3 had no such oscillation, 
but 
increased with drop size. 


possessed a degree of oblateness which 


Drops for all three 








cm/sec 


| 
: 
> 











curves fell with their major axis normal to the 
direction of motion. 


Effects of oblateness on fall velocity 


Curves 2 and 3 fall further below the predicted 
curve for rigid spheres as drop size increases, 
presumably due to increasing form drag with the 
observed increase in oblateness. 


Effects of internal circulation on fall velocity 

The weak circulation within drops of curve 2 
evidently reduces the total drag sufficiently to 
cause all these drops to fall some 2-4% faster 
than the stagnant drops of curve 3. 

The unshaded points of curve 1 correspond to 
drops of volume c.a. 0-017cce; they showed 
vigorous internal circulation but no detectable 
They fall about 40% faster than 
predicted for a rigid sphere of the same density 
and size, and about 50% faster than the drop of 
the same density and size in curve 3, this latter 
drop being completely stagnant. 


oscillation. 


Effects of oscillation on fall velocity 

Since all 4 curves must pass through the origin, 
then a maximum fall velocity is reached for drops 
of curve 1, below a drop size of 0-02 cc. Above 
this drop size oscillation of the drop occurs, 
reducing the fall velocity. The extent of this 
oscillation increases with drop size. If we assume 
that, in the absence of oscillation, curve 1 would 
have flattened off at larger drop sizes in a similar 
manner to curve 2 or 8, then curve 1 shows that, 
for a 0-08 ce drop, the oscillation occurring in this 
system reduced the fall velocity by at least 12%. 
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On the steady state fractionation of multicomponent and complex mixtures 
in an ideal cascade 


Part IV. The calculation of the minimum reflux ratio 
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Summary—In this paper formulae for calculating the minimum reflux in rectification are 
summarized and developed for complex mixtures which are characterised by continuous or 


semi-continuous distribution functions. 


Two numerical examples are considered, one for 


continuous and one for semi-continuous mixtures. 


Résumé—Dans cet article, les auteurs résument et développent des formules permettant le 
calcul du reflux minimum dans la rectification de mélanges complexes caractérisés par des fonctions 


de distribution continues ou semi-continues. 


Ils examinent deux examples numériques : un 


pour les mélanges continus et un pour les mélanges semi-continus. 


INTRODUCTION 


In a preceding paper the present authors presented 
a solution, in closed form, of the equations for 
the separation of complex mixtures, that is, 
mixtures with an indefinite number of components, 
in an ideal column. In this paper will be 
summarized the solution to the minimum reflux 
problem for complex mixtures, together with two 
illustrative examples. 

First of all, it is to be noted that at any given 
reflux ratio the separation of a feed of known 
composition and heat content is completely 
specified if two and only two points of the 
distribution function of the overhead product 
are fixed a priori, provided, of course, that the 
feed is introduced at the optimum location in the 
column. The above is but another formulation 
of the well-known theorem that, for any multi- 
component or complex mixture, the number 
of degrees of freedom in a distillation column 
is equal to 2. Secondly, it will be recalled that 
for a specified separation of a given feed, as the 
reflux ratio approaches its minimum permissible 
value, both N and S, the number of plates in 


the enriching and the stripping sections of the 
column, respectively, increase beyond limit. 
Finally, use will be made of the following theorem, 
which was rigorously proved, and which forms 
the basis for this discussion of the minimum 


reflux problem. 
THEOREM: 


At minimum reflux, and for a specified separation 
of a feed of known composition and heat content, 
there exist two values of p, namely p = h and 
p = l, such that : 


(a) a(p) > 0, if ap(p) >0, for p<h 


N 
(b) lim tolP) (") exists for p >h 


N->@ Ty(h) \h 


(c) x* (p) > 0, if zp(p) > 90, for p>l 
a*(p) (: 

P 
where 2o(p), @p(p) and a*(p) are, respectively, 
the distribution functions of the overhead 


product, the feed and the bottom product. 
Following a well-established tradition, the point 


(d) lim 


So as (i) 


s 
) exists for p< l 


* Department of Chemistry and Chemical Engineering, University of California, Berkeley, California. 
** On Leave, Cambridge University, Cambridge, England. 
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h will be referred to as the heavy key component 
and | as the light key component. 

The minimum reflux problem, then, under the 
assumptions of ideal plates, constant molal 
overflow and constant relative volatilities through- 
out the column, can be solved exactly for con- 
tinuous mixtures as follows : 


First it is supposed, for simplicity, that the 
distribution function z,{ p) is everywhere positive 
on the positive p-axis. Then the following 
functions are defined : 


(a) wed) = —q— P[ PAP) ap forl <A <h (1) 
0 + 


where the symbol P in front of an improper 
integral denotes its Cauchy principal value and 
q represents the thermal condition of the feed. 

_) 7 AAA) 


4 
(b) (A) = - tan i 


(c) Bip) = 


for! SA sh (2) 


h 
o(p) afr | Aaa —¥ [e(r) |" enh (3) 
P ; ; 


for! Sp sh, where { (2i) is the Riemann zeta 


function [1]. However, if p(p) is close to unity, 
then it is better to rearrange eq. (3) into the more 
convenient form 


B(p) = op) [1 — (p(p))*] 
exp ‘" a ae (Pi (a2) - | (4) 
—_ 1 
i 


for! Sp sh. 


(a) Bip) = exp f 


r 


[2 | forp<l (5) 
A—p 
i 


. i 
Az pA) 
Dpzdp) = Bp) {eo + | wok — al (6) 
J BO) (A—p) 
for! Sp Sh, and 
Az (A) 


8= } Bay) 


dd (7) 


while 
@ t h 
D = | Dad) dp =| =X) dp +[Dalp) dp (8) 
i 


0 0 


Equations (6), (7) and (8) enable one to calculate 
not only the minimum reflux ratio but the dis- 
tribution function of the overhead product as 
well. Moreover, the relations 


i h 
Rp — | 22H?) 4 + [PpaePa 9 
vw P Fees (9) 


and 


7 


Dpzp) = ~ tan mp(p) {RD - 


pf Ded) 4, .. SP) ran ols) 


— 4 (10) 


forl sp sh 


can be used as a check in any numerical 
calculation. 

It is clear from the above formulae that 
the minimum reflux problem, for any desired 
separation of the feed, can be easily solved 
if the two points A and I, that is, the two key 
components, can be spotted in advance. Since 
this cannot be accomplished readily in most 
cases, it appears preferable to assume values for 
h and / until the calculated separation of the 
feed becomes essentially identical with that 
specified in the problem. 

This brief presentation of the analytic method 
for solving the minimum reflux problem has 
been restricted to continuous mixtures for which 
the distribution function was positive along 
the positive p-axis. 


ILLUSTRATIVE EXAMPLES 


(a) Continuous Miztures. 


Let us suppose that the distribution function 
of the feed, z,{ p) is given by 


zp) = a(P — 1}%e--"" forall p >1 


and z,(p) = 0 for p< 1. The vapour-to-liquid 
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Table II. Numerical values for w A), tan mp(A) and p(A). 


tan wm p (A) p (A) wp (A) | tan wp (A) 





0-5318 0-1556 1-5878 | — 3-2026 
0-6401 0-1812 19376 | — 2-5329 
0-7626 0-2074 2-2562 — 2-0821 
0-9034 0-2339 2-5406 — 1°7560 
1-0686 0-2606 2-7888 — 1-5072 
1-2674 0-2874 3-0000 — 13102 
15143 0-3142 38-1741 — 11493 
1-8332 0-3410 33119 | — 10149 
2-2654 0-3677 3-4151 — 0-9003 
29015 | 060-3944 38-4856 — 0-8012 
38-9327 0-4207 35261 ~ 0-7144 
5-9358 0-4469 35394 0-6376 
11-6400 0-4727 3-5286 0-5690 
181-4916 0-4983 3-4970 0-5075 
— 13-5750 0-5234 — 34477 | 0-4520 
— 65588 0-5482 - 39-3837 | 0-4017 
— 43154 0-5725 




















Table III. Numerical values of p(A) at shorter ,(B) 
intervals of A x*(B) 


r 
= 20 and - 





aiso spec * é 
oa) oa) a) It is also specified that 
— h = 3-05 and / = 1-45. 

0-1556 0-4469 0-7075 
0-1683 0-4598 0-7178 The above is a somewhat simplified problem, 
0-1812 0-4727 07280 since one usually does not specify h and 1, but 
0-1943 0-4856 0-7378 
0-2074 0-4983 0-7476 
0-2206 0-5109 | 0-7573 Table IV. Numerical values for B(A), 1-45 
0-2339 0-5234 0-7667 <A<805 
0-2472 0-5359 0-7759 
0-2606 0-5482 0-7850 
0-2740 0-5604 0-7939 B (A) B (A) 
0-2874 0-5725 0-8026 = 
0-3008 0-5845 0-8110 ‘ 0-72012 0-59574 
0-3142 0-5964 0-8193 | 0°76360 “ 0°53131 
0-3276 0-608 1 0-8274 081799 046740 
0-3410 0-6197 ’ 0-8353 0-86952 0-40611 
0-3543 0-6313 0-8430 d 0-91224 0-34774 
0-3677 0-6426 0-8505 | 094416 0-29324 
0-3811 0-6538 0-8578 ’ 096436 ; 0-24303 
0-8944 . 0-6648 0-8649 ‘ 0-97220 0-19808 
0-4075 0-6758 0-8717 096659 0-15772 
0-4207 . 0-6865 | 0-8784 0-94956 0-12238 
0-4339 0-6970 | 0-92096 . 0-09186 
| | 0-88237 0-066038 
0-04400 
0-02664 
































ratio in the feed is equal to 1. It is desired to 
determine the minimum reflux ratio and the two 
points p = f and p = f* such that 
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B and §* instead; however, it will suffice as an 
example. Some of the results are presented in 
Tables I-V. Table I will be found in the previous 
paper, Part 3. The method of calculating 
Cauchy principle values of integrals was also 
discussed in Part 3. 

It can be shown now from eq. (3) that in any 
small interval AA, the function 


Arg(A) 


BQ) (1-45 — A) 


Table V. Numerical values for B (A) and 
Arp(A) 


— =A") __ for 1.00 < A < 1-45. 
Bay 0-45 — ay” 








0-52696 
0-51287 
0-49889 
0-48503 
0-43077 
0-33019 
0-24132 
0-16563 
0-10408 
0-05710 
0-02458 
0-00591 


| 467410 
432640 
| 410071 
| 854360 
3-02391 
2-72976 
2-52797 
2-37669 
2-25696 
2-15937 
207773 


4°73120 
2°43120 
1-09193 
0-58935 
0-82760 
O-17517 
0-08433 
0-03253 
0-00711 











can be represented as K/(1-45 — A)* where both 
K and « vary slowly with A. It is found then 
that for 

1-44875 SA 51-45, a = 0.8444 = 1 — p(1-45) 
1-4375 SA 5 1-44875, a = 0.86656 

1425 SAS 1-4375, « = 0-89174 

1.40 SAS1-425, «= 0-95313 


Therefore, from eq. (7) 


1-45 
“4 Aw,(A) 
-| B(A) (1-45 — A) 
The function p Dap) can next be calculated 
from eq. (6) and must be identical with that 
given by eq. (10), namely 


AY) ton eof) 


dd = 1.04971 


[A Dro A) Jeate. = 


This comparison is made in Table VI. 
Finally, from eq. (9) 


(RD)... = 1.04895 
Also 


@ 


D = | D2,(p) dp = 0-50991 


ri 
and therefore 


Rain, = 20586 





A Daf A) 











0-60968 
0-54445 
0-47950 
0-41706 
0-35745 
0-30169 
0-25022 
0-20410 
0-16262 
0-12628 
0-09484 
0-06820 
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It is seen then from Table VI that, on the 
whole, the agreement between A Dz,(A) and 
[A Dao(A) aes a8 well as between RD and 
(RD)..;., is very satisfactory. It was hoped 
moreover that by means of eq. (10) one could 
have obtained even more accurate values for 
A Da,A) by applying an iterative method. 
Various attempts to discover the appropriate 
iteration were not successful, however, and the 
Seidel method, which is the most commonly 
used in such cases, did seem to diverge, in spite 
of the close agreement between A Da,(A) and 
[A Dao(A) |i... It appears therefore that the 
best way to improve the accuracy of the results 
would be to increase the accuracy of the previous 
calculations. 


It is necessary to obtain the two components 
B and 8*. This can be done with the aid of 
Table VI, and it is estimated that 


B ~ 1-495 
and 
p* ~ 2.9524 


In this problem then the two “key com- 
ponents "’ were specified a priori and then £ 
and £* were determined ; it is much more common, 
however, in practice to specify 8 and §*, and then 
to calculate the minimum reflux. The procedure 
adopted here however appears quite straight- 
forward, and therefore it should be used even 
when A and / are not known. In other words, 
when f and §* are specified, one assumes different 
values for h and / and evaluates the corresponding 
Ri,» 8 and B*; the desired answer can then be 
obtained by a proper interpolation. 


(b) Semi-continuous Miztures. 

It is supposed that the distribution function of 
the feed, z,(p), is now given by 

zp) = a (p —1)*e-*-"" for all p>1 
andz,p) = 0 for p < 1. It is also supposed that 
one component is present, at p = 2 50, in finite 
amounts, and that its mole fraction in the feed, 
a,li), is equal to 0-25. The vapour to liquid 
ratio in the feed is again taken equal to 1. It is 


desired to determine the minimum reflux ratio, 
and the two points p = § and p = §*, such that 


wo(B) _ x5(8*) = 20 
74(6) z.(B*) 
It is also specified that 
h = 8-05 and | = 1-45 


This problem then is quite similar to that worked 
out earlier, except for the presence of the dis- 
tinguishable component at p = 2-50. The func- 
tion p(A) must now be determined from 


AxgA) 2 cot mp(A) = — 0-50— 
ra Pg xxhq) —_ 
P| f-2 dq = w,(A) 


10 


0-6250 
A — 2-50 


_y 7 A@,(A) 
wp A) 
Numerical results are given in Tables VII-X. 
So far, equations (1)-(10) have been used with 
only one modification. However, the following 
changes are necessary, now that one distinguish- 
able component is present in the mixture : 


(a) 


p{A) ax - tan 


Instead of eq. (8) 


= [ Dedp)dp + Dag) (a0 
oO 


(b) Instead of eqs. (9) and (10), respectively, 
8.08 


Rp - { p D2fp) , 

8-05 — p 
0 

and 


2.50 D z,{i) 


12 
8-05 — 2-50 (12) 


D pap) = ~ tan vom - 


~P (Pea 


2-50 D xi) 
2-50 — p Pp 


a = Se) (18) 


tan mop) for! sp sh 
where 2,(i) is the mole fraction of the distin- 
guishable component in the overhead product 
and is given by 
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Table VII. Numerical values for w,{A), tan = p (A) and p(A). 





| 
w pf A) tan pf A) pt A) w p{ A) | tan mp A) 


38238 0-4446 0-1332 2-4122 2-1081 
39014 0-5307 | 1553 34512 14220 
3-9319 0-6248 0-1777 35-9104 0-7948 
3-9132 0-7281 0-2003 13-9594 0-3196 
38453 0-8424 2228 
3-7298 0-9700 2452 -~19833 | 0-1982 

3-5706 11137 2671 ~ 11-6741 0-3125 

3-3731 1-2768 2885 9-3721 0-0342 0-8866 
31460 14626 3091 2 - 7-7484 - 03968 | 0-8798 
2-8951 1 -6766 3288 75 6-9856 0-3998 0-8789 
2-6352 1-9203 3472 64705 0-3892 «= |:«O-8818 
2-3788 2-1929 3638 6-0870 0-3707 0-8870 
2-1423 2-4831 0-3782 57786 0-3475 0-8936 
19463 2-7603 0-3894 5-5155 0-3218 0-9009 
1-8186 2-9538 0-3961 5-2810 0-2951 0-9087 
1-7997 2-9569 0-3962 5-0655 - 02683 0-9166 
1-9559 2-6712 0-3860 








VIII. Numerical values at shorter Table 1X. Numerical values for 
intervals of A. <A< 3-05. 





p(A) p{A) p( A) | BOA) 
0-1332 . 0-3638 5! 0-9377 “! 0-50599 : 0-14282 
0-1442 . 0-3714 57! 0-9185 ‘55 | 0-56193 0-60266 0-11341 
0-1553 | 2 0-3782 4 0-9036 0-58356 0-56387 7 0-08774 
01665 | 2075 o3s63 | 2-62: 0-8938 6: 0-62204 0-52136 ‘80 | 0-06568 
0-1777 : 0-3894 6: 0-8866 0-65420 | O-47614 . 0-04701 
0-1890 : 0-3936 67! O-8824 ‘7 0-67825 | O-42807 0-03144 
0-2003 SF 0-3961 0-8798 0-69353 |‘ 0-38081 . 0-01868 
0-2116 47! 0-3974 0-8789 8: 0-699044 | 240 | 0-33607 ;, 0-00842 
0-2228 2 0-3962 0-8789 0-69622 020813 
0-2340 . 0-3932 0-8801 , | (68417 2°: 021817 
0-2452 | 2-2! 0-3860 0-8818 . 0-66383 | 2- 0-17663 
0-2562 2-27: 0-3759 O-8842 

0-2671 4 0-3590 0-8870 
0-2779 “f 0-3366 87! 0-8902 re 
0-2885 s 03049 2-900 0-8936 So, it is found from eq. (7) that 
0-2989 0-2624 | 2-925 0-8972 1-48 

0-3001 ; 0-2138 2-050 0-9009 R | q 7 lq) 
0-3191 “3 0-15.40 2-975 0-9048 Biq) (1-45 — q) 
0-3288 . 0-0985 | 3-000 0-9087 
“Be » . of q 3° 25 ’ « 
pom ae eo oe Next, the function p Dag(p) can be calculated 
0-3558 | 2-525 | 0-9653 from eq. (6), and must, moreover, be identical 


tS t bt 


~ 


te 











~ 


dg = 1-10911 














with that obtained from eq. (13), namely 


zi) = Hi) 2 = 0.25 Bs [A Dr) Jeate. = = tan mp(A) 
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Table X. Numerical values for B(A) and 


Arp(A) 
Rtas — Taal 1.00 <A < 1-45. 








Aw p(A)/B (A) (1-45 — A) 


Aen) | BO) 





8-21424 
4°34300 
2-92893 
2-17640 
0-96158 
0-510738 
0°28043 
0-14850 
0-07091 
0-02716 
0-00590 


3-74617 
3-35040 
8-12677 
296894 
2-57535 
2-36251 
221483 
2-10266 
2-01340 
1-93962 
1-87789 


038465 
036377 
0-84343 
0-32308 
024764 
0-18099 
0-12422 
0-07806 
0-04283 
0-01844 
0-00443 














and [A Dz,(A)]..:., 28 well as between RD and 
(RD)... is very satisfactory. As was said 
earlier, in connection with the first problem, 
the agreement between the two functions could 
not be improved by any interative method, 
and again it was found that the Seidel method 
diverged. So it appears that more accurate 
results could be obtained only by increasing the 
accuracy of the calculations. 


Finally, the two components 8 and £* can be 
estimated from Table XI. It is found that 
B ~ 1-406 
and p* ~ 2-949 


It is hoped then that the example just presented 


Table XI. 





A Darl) [A Dizol A) Jeatc. 


A Dag A) [A Dato A) Jeate. 





0-47925 0-47867 
0-54048 
0-59566 
0-64453 
0-68483 
0-71534 
0°-73565 
0-74487 
0°74459 
0-73394 
0-71895 
0-68587 
0-65083 
0-60995 
0-56479 
0-51647 


0-59480 


068471 


0-73597 


0°74637 


0-71624 


065441 


0-56967 











046554 0-47127 
0-41399 
0-36570 
0-32469 
0-27470 estim. 
0-23797 
0-19277 
0-15598 
0-12393 
0-09593 
0-07304 
0-05144 
0-03443 
0-02046 
0-00922 


0-37492 











This comparison is made in Table XI. 
We also have that 


To (2-5) 
x p( 2-5) 


oli) = 0-25 = 0.068436 
and that 
(RD),.,.. = 1-11888 
Finally, D = 0-4588, and therefore, 
Rin. = 2-420 


It is seen then from Table XI that again on 
the whole, the agreement between A Dzo{A) 


Table XII. Some numerical values 
of the Riemann zeta function. 





{ (2%) 





1-644934 
1-082323 
1017343 
1-004077 
1-000995 
1-000246 
1-000061 
1-000015 
1-000004 





oe ouek @ Ww 
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will adequately illustrate the semi-continuous at least, be essentially identical to the one just 
problem; because the problem with more discussed. Table XII gives some numerical 
distinguishable components would, in principle values for the Riemann zeta function. 
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Dimensional systems and systems of units in physics with special 
reference to chemical engineering 


Part II. 


Practical rules for authors and readers and the need for their observance 


A. KLINKENBERG 
N.V. De Bataafsche Petroleum Maatschappij (Royal Dutch/Shell Group), The Hague 


( Received 23 January 1955) 


1. INTRODUCTION AND SUMMARY 


In Part I the principles were discussed of con- 
structing systems of units. The results of this 
study will now be applied to formulate a number 
of practical rules, the application of which should 
reduce the amount of misunderstanding in the 
field. 

Most of these rules are intended for authors 
(see § 2). The main line of thought therein refers 
to the case where an equation in physics or 
engineering is derived from first principles and 
then numerically applied (rules 1-9). 

Separate statements are given with reference to 


empirical equations (rule 10), 

dimensional analysis (rule 11), 

dimensionless groups (rule 12), 

equations that should be applicable in more 
than one dimensional system (rule 13). 


In § 3 a few recommendations are given for the 
use of the readers. 

In § 4 a number of literature quotations are 
analysed. This discussion shows that many 
authors put problems before the readers, which, 
it is suggested, they could have avoided. 

Finally, a summary is presented of advantages 
and disadvantages of various systems, either 
derived by a priori considerations or from the 
analysis of the literature. 

The weighing of these arguments is of course a 
subjective procedure. It is considered that 


following systems are preferable : 
Metric: the practical system (m, kg (mass), 
sec., °C, practical electrical units), 

the Ib (mass), lb (force), ft., sec., BTU, 


°F system. 


British : 


RULES FOR THE AUTHOR 


1. The advice is given to consider letter symbols to 
stand for physical magnitudes and nut for their 
numerical values only. 


2. PRACTICAL 


This means that statements like : “ the accelera- 
tion by gravity is g m/sec.” ” are not allowed. 


2. Decide on a coherent set of fundamental quantities 
and dimensional constants and use it consistently. 


The sets of quantities discussed in Part I are : 


In mechanics : 

(a) MLT* : dynamic systems. 

(6) FLT: static systems. 

(c) FMLT: English engineering system. 


In mechanics + heat : 

(d) MLTQ®, e.g.: cgs system : g,,, cm, sec., cal., 
°C. Brit. dyn. system : Ib.,, ft., sec., BTU, 
°F. 

(e) MLT®, e.g.: pract. system: kg,,, m, sec., °C. 

(f) FLTQ®, e.g.: kg,, m, sec., keal., °C. 
lb.,, ft., sec., BTU, °F. 

(g) FMLTQ®, e.g.: lb.,, lb.,,, ft., sec., BTU, °F. 


In mechanics + heat + electricity : 

(h) MLTQ®, the cgs system: g,,, cm, sec., cal., 
°C, 

(i) MLTOI, the practical system : kg,,, m, sec., 
°C, ampére. 


The above list represents a choice of the most 
important combinations of wide applicability. 
Others, such as practical electrical units + 
mechanical cgs units, are also met. 


* M = Mass, L =Length, T= Time, F = Force, 
Q = Heat, @ = Temperature, I = Electric current. 
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Excluding applications of Newton’s First Law 


(mutual attraction of masses), the dimensional 


constants to be used with the above sets of 


quantities are : 


In mechanics : 

(a) none. 

(b) none. 

[ in Ib., Ib... ft. sec. system : 
_ og Ib, ft. 

} geo Ib., sec.” | 

in Ib., Ib.,, ft. hr. system : 


ML 
(c) & = fe 





In mechanics + heat : 

27-2 
(d) J= = T 
Q 


= mech. equivalent of heat 


in egs system : 


J = 419 x 107 “8 


cal 
in Brit. dyn. system : 
I= %1 ft. ae 
(e) none. 


in lb., ft. sec. BTU. “F system : 


J = 778 a 
BTU 
fe as under (c) 

FT?|’ 


in Ib. Ib.,, ft. see. BTU. °F — 
Ib., ft. 


778 —! 


BTU 


system : J = 


In mechanics + heat + electricity : 
(h) one dimensional constant in some formulae in 
electricity and magnetism. 
; = 4n V sec. 
i = permeability of vacuum = —. — 
( ) Ho pe . 10° A m 
€y = dielectric constant of vacuum - 
10’ A sec. 9-086 x 19°? A sec. 
~ 4n(8 x 10°? Vm ; Vm 
Examples A 11 and A 12 in § 4 show difficulties 
that arise through unsystematic use of dimensional 


systems. 


8. Decide on factors firing numerical constants. 
Some of the alternatives have been mentioned in 
Part I, e.g. : 
rationalization or non-rationalization in the 
field of electricity and magnetism ; 
various ways of defining a Reynolds group 
(radius vs. diameter of pipe ; particle size vs. 
inverse of specific particle surface) ; 
various ways of defining the friction factor for 


flow through a pipe ( the American : f 


shear stress at wall divided by }pv*, the 
European : A = 4f = pressure drop divided 


by dpv*. ‘): 


Deviations from accepted custom in 
respect cannot be traced by dimensional checks. 

It is impossible to give precise instructions. 
Rule 3 only serves as a reminder that in deriving 
physical relationships a number of conventions 
are used, often unconsciously. 


this 


4. Prepare a list of quantities and dimensional con- 
stants giving their symbols, names and dimen- 
stonal formulae, and keep this list complete. 


. Write down equations in letter form, ensuring 
dimensional homogeneity by using appropriate 
dimensional constants. Carry out the necessary 
mathematical transformations and check dimen- 
stonal homogeneity of the final result. 


Remark—The final check is in principle super- 
fluous but it is very useful to help in detecting 
mistakes. 

6. With equations in letter form, do not prescribe the 
units to be used when this constitutes an unneces- 
sary limitation. 

. Choose a unit for each fundamental quantity. 
Substitute the units in the dimensional formulae 
for the derived quantities in order to obtain the 
units of these derived quantities. 

. In order to prepare the equation for numerical use, 
replace constants (dimensional constants, phy- 
sical quantities that happen to be constant in 
the given problem) by their numerical values. 
From this point onward, mention each time the 
units to be used. 
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It is convenient to mention the unit of the 
numerical factor in order that users of the equa- 
tion may easily convert it into their own system. 

THORNE Watsnaw in their booklet 
“ Engineering Units and the Stroud Convention ” 
(1954) give prescriptions how to derive numerical 


and 


equations from physical equations. 

Their physical equations, which cover the fields 
of mechanics and heat are always written without 
the dimensional constants, g, and J. This is how 
they come to the statement (p. 2) that physical 
equations are completely free from any restric- 
tions as regards particular units or systems of 
units (see, however, Part I, page 132). 

In a physical 
equation THorne and Wa.tsnaw then replace the 
letter symbols by products of numerics and units. 
No prescription to use a coherent system of units 


dimensionally homogeneous 


is given so that the resultant equation is “ uni- 
tarily inhomogeneous.” By inspection the so- 


‘unity brackets ” 


called * are then chosen and 
added as extra factors (each being unity) to 
obtain unitary homogeneity. These factors 


cover : 


the conversion of size of unit : 
5°C 3600 sec. 
e.g.: 3 
oe 1 hr. 
the replacement of derived units by fundamental 


units : 
gE sec.” E sec.” 
e.g. : — I; . 
gem slug ft. 


the change in fundamental unit : 


i Ib.,, ft.] & ft. Ib.¢] F= hr. 
 Tb.psec.? BTU 6080 ft. 


, 


This subdivision of “ unity brackets ” in three 
groups has been mentioned here for the purpose 
of seeing the relations with the present article. 
It is, however, artificial and we may just as well 
say that by not having systems of units, THorNe 
and WaLsHaAw leave the equations in physics full 
of “ dimensional constants,” which they do not 
write in equations in letter form and must intro- 
duce later. 


9. Give the list sub 4 at the end of the article. If 


desired the dimensional formulae can be replaced 
by the units. 
Remarks- 
alternatives. 


1. Obviously, such a list cannot allow 


2. In expressing the units it is not necessary to 
use fundamental units only. It may be more con- 
venient to use some very important derived units, 
such as: dyne, erg, newton, joule, or, volt. They 
have precisely been given their names for that 
that purpose. 


10. Empirical equations that are not dimensionally 
homogeneous, should never be given without 
stating the units of the various letter factors. 
Preferably the unit of any dimensional numer- 
ical factor should be mentioned in order to 
facilitate conversion. 


Dimensional analysis need not be carried out 
with the fundamental units the author uses gen- 
erally. There is an advantage in using the 
maximum number of fundamental units com- 
patible with the problem. No unnecessary 
dimensional constants must be included in 
the set of quantities to which dimensional 
analysis is applied. 
Remarks—1. In heat transfer without conversion 
of heat in mechanical energy or vice versa, heat 
will be given a separate dimension —even if 
practical units are used elsewhere — without, how- 
ever, introducing the mechanical equivalent of 


heat. 


2. A statement by GuGcENHEI™M [1] is highly 
appropriate : 


“. . . the number of fundamental quantities 
having independent dimensions is, to some extent, 
a matter of choice. But, if in the same problem 
or set of problems two authors make a different 
choice, the one choosing the greater number is 
likely to be the more competent physicist.” 


12. Use dimensionless groups to the utmust. 


All complete equations can be written as relations 
between dimensionless groups. 

There is an advantage, in giving data and in 
presenting relationships in analytical or graphical 
form to use dimensionless groups to the greatest 
possible extent, since : 
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. The numerical values of dimensionless groups 
are independent of the choice of the system of 
dimensions and the further choice of the sizes 
of the units ; 


. The number of variables is reduced thereby, so 
that equations assume their simplest form and 
the amount of experimentation to derive an 
empirical relationship is a minimum. 

It is advantageous to introduce dimensionless 
groups at an early stage of the mathematical 
transformation, for instance to write already the 
governing differential equation and boundary 
conditions in dimensionless form. 

In calculating numerical values of dimensionless 
groups coherent units must be used. Thus, if each 
symbol for a physical magnitude is replaced by the 
product of a numerical value and a unit, expressed 
in fundamental units, all units cancel. 

Examples B 5 and 7-13 signal cases where the 
units as prescribed are inconsistent and yet the 
correct result is given. Thus, the author did not 
follow his own prescription. 

In Examples B1 and 6 the units are not 
coherent and the result deviates from the generally 


accepted one. 


13. Multi-purpose formulae. 


If formulae are to be ready for use with more than 

one system of units, following rules apply : 

(a) If all systems of units, in the field covered by 
the calculations, require the same numerical 
factors and the 
equations in letter form are valid for all of 
them. 


dimensional constants, 


If numerical factors and /or dimensional con- 
stants differ, the equations should be provided 
with factors and constants for all cases, with 
precise instructions where to drop those. 
The survey by Kronic [2] of equations for 
rationalized practical units and for various types 
of egs electrical units is given in this manner. 
The combination of a dynamical system with 
the English Engineering system would require the 
introduction of g, with the instruction to drop it 
for the dynamical systems. Failure to do so 
cannot be corrected by conversion of units (vide 
example A 8). 


8. Practica RULES FoR THE READER 
1. Check on observance of rules by author. 


2. If desired, convert equations in letter form to the 
own preferred system by following methods : 

(a) If the set of dimensional constants in the 

reader's preferred system is comprised in the set 


used by author : 


Omit any superfluous dimensional constants. 


Example : 
Author: g.F = ma (Engl. Engineering). 
Reader: F = ma (static and dynamic systems). 


(b) If the reader’s set contains new dimensional 
constants : 
Find by a dimensional check which dimensional con- 
stants are to be included. 


Example : 
Author: pV = RT (practical system). 

p = (ML T*). 

V = [L* Mol"). 

R = (ML? T? Mol" 6} 

T = (9). 
J RT (cegs system + cal + °K). 

p = (ML" T*). 

V = [L* Mot"). 

R = [QMol" 6"}. 

T = (9). 

J = (ML? T? Q"). 

Remark—Caution is necessary in the fields of 

electricity and magnetism in converting cgs into 
practical units and vice versa. Equations may 
differ by a numerical factor 47, due to the fact 
that the practical system is rationalized and the 
egs system is as a rule not. In cases of doubt, it 
is recommended to consult tables. 


Reader: pV 


so that : 


3. Conversion of units of physical quantities and 
constants in numerical equations. 


(a) Conversion within the same dimensional system. 


If each quantity is considered to be the product 
of a numerical factor and a unit the conversion 
is easily carried out with the aid of the conversion 
factors for the fundamental units. 

One method is, simply to replace each funda- 
mental unit from one system by its equivalent in 
the other. Thus, for example, 
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100 Ibs. /cu.ft. = 
= 100 x 0-454 kg/(0-3048 m)* = 1602 kg/m* 


If a ready conversion factor is used, it may be 
applied in the wrong direction unless one uses 
the principle adopted by van Disck [3] of heading 
tables in the following manner : 





| Er- | Ex- 


> to multiply with + 
d\~ te 
Quantity _ + to divide by «- pressed 


| 
force | 





dyne | newton 





(b) Conversion to another dimensional system. 


It is probably safer not to attempt to convert to 
units of other dimensional systems but to change 
the dimensional system already in the equation 
in letter symbols, as mentioned in para. 2. 

“xample A 8 indicates the sort of difficulties 
one might otherwise encounter. 


4. EXAMPLES 


In most cases the authors make no numerical 
mistakes. Yet their prescriptions, if duly carried 
out, would often lead one to do so. 

It is the aim of the following discussion to prove 
that insufficient attention in this respect places 
an extra burden on the reader. 

The examples are grouped under the following 
headings : 

A. Unsystematic use of dimensional formulae 
and dimensional constants 
statements thereon. 


and incorrect 
Inconsistencies in the substitution of numeri- 
cal values in equations. 
Regarding subject B, the substitution of 
numerical values in equations, it is noted that the 
author may : 


. give an incomplete description of the units to 
be used in his equations ; 

. give an overcomplete description, i.e. leave 
alternatives to choose from ; 

. describe an inconsistent set of units. 

In this respect it is most dangerous to omit a 

dimensional constant from a coherent set, when 
the incoherence easily escapes attention ; 


4. describe no units at all. 

This is in many cases acceptable, viz. when 
no numerical application is given, when it is 
clear what dimensional system has been used 
and it is left to the reader to choose the sizes 
of the fundamental units. 


. Unsystematic use of dimensional formulae and 
dimensional constants, incorrect statements on 
dimensional systems. 


Al. Author Al* under a_ heading 
system’ mentions erroneously : 


“ 


cgs 


conversion factor g, = 980-6 8m Cm 
&r 


unit for the gas constant R : &eem 
mol °C 

In a column “ English” analogous definitions 
are given together with a statement that the unit 
of force is the poundal. However, lb.,and poundal 
are incompatible. 

Elsewhere g. is incorrectly given as 82-17 
(Ib... /Ib.,)/(ft. /see.). 


A2. Authors A 2 in their notation use F, M, 
L, T, Qand 6 as fundamental quantities and give 
corresponding dimensional formulae, such as 
pressure drop [FL™] and density = [ML™*]. This 
is the English Engineering system, new version. 

The text, however, equates a pressure drop to 
a density times a length, which means that F and 
M have been equated, as described on page 136, 
Part I (English Engineering system, old version.) 


A3 and 4. At one place the author of A3 
presents an admirable treatise on the use of the 
English Engineering system (new version) with 
the dimensional constant g.. At other places 
there are some misstatements regarding the unit 
in which this constant is expressed, viz.: 


g. = 418 x 10° lb force x 
ft. /(pounds fluid) (hr.) (hr.) 


A 3, elsewhere, 


Ib matter ft. 
hr.? 


The authors of A5 make no distinct 


A 4, same author, g, = 4:17 x 10° 


A 5. 


* The list of authors has been submitted to the Editor 
of this Journal but will not be published. 
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choice between the three systems mentioned in 
§ 2, rule 2, a, b, and c, having 


fundamental units derived units 
(a) 'b.,, ft., see. 
(b) Ib, ft., see. 


(c) Ib... Ib.,, ft., see. 


poundal 
slug 


They define four fundamental dimensions 
(system (c)) but use the units sub. (6) for them, 
one of which is a derived one. Newton's Law is 
then used to derive a relation between the funda- 
mental! dimensions (logically impossible!) in order 
to decrease their number to three (system (a) or 


(b)). 


The notation is unsystematic, e.g. 


ft.lb., sec. 
slug, cu.ft. or 
Ib. cu.ft. 
viscosity Ib... /ft.sec. 
Ib. 

slug a aan (b) 


power occurs in system (b) and (c) 


density (6) or 

(a) and (c) 
(a) and (c) 
force (6) and (c) 


mass 


A 6, 


a porous bed measuring pressure in g cm* ; 


Author A 6 defines a permeability A, of 
vis- 
cosity in poises, length in em and linear velocity 
in cm/sec. Thus the quantity defined is the pro- 
duct of the permeability in cm? (usual definition) 
and the dimensional constant g in cm ‘sec.” 

As a result of the above, an equation is obtained 
relating permeability to porosity (dimensionless) 
and specific surface S, (em™') of following form : 


8 
K, a af . 2 

5S) (1 — €) 
where g makes a mysterious appearance. 

The List of Symbols says that all units are 
consistent cgs units except where otherwise stated. 
It omits to give a unit for K,, and this is definitely 
not a cgs unit. 


A7. Authors A7 give author’s A 6 equation 
in the following form : 


S, ere 
eae 


Commenting thereon, in this equation the 


factor 14 is a dimensional constant ( nT s} Again 


the dimension of K, is omitted in the Notation. 
The equation in which K, is defined is not of 
much avail, since the prescription is given to insert 
the viscosity in centipoises instead of poises. 

A8. Authors A8 use formulae Valid in 
dynamical systems, i.e. not including g.. They 
wanted to make the result applicable with the 
English Engineering system (using hours) as well, 
by indicating conversion factors in the list of 
symbols. 

The pressure drop in N'm?* then has to be 
replaced by g, times the pressure drop in |by/sq.ft. 
To say Ap = pressure drop N/m? (= 8-7 = 10° 
lb. force /ft.*) is incorrect. 

A9%. Author A 9 gives in one article formulae 
valid for the metric equivalent of the English 
Engineering system (old version), involving a 
dimensional constant g (in em.sec.~*), and formulae 
valid for the cgs system. The Table of Symbols 
says that pressure is expressed in cm water gauge 
or gem™, which is incompatible with the cgs 
formulae. 


A 10. Authors A 10 have expressed the pressure 
difference over a pipe section, Ap, in static units 
(Ib ‘sq.ft.) and the shearing stress F at the wall 
of the pipe in dynamic units (Ib ft.hr.*) in accord- 
ance with their list of symbols. When using the 
English Engineering system a choice must be 
made, as mentioned in Part I (page 136) whether 
quantities are expressed in the one way or the 
other. It is very unfortunate that in defining the 
units for two so closely related quantities as a 
normal and a tangential force per unit area, 
different routes have been used. 

Consequently the equation relating Ap, F and 
the pipe dimensions (d and 1) contains g as a 
dimensional constant, viz.: 

nd F 
<< woe 
with g = 418 = 10° ft. /hr.? 

However, later on, again for pressure drop along 
a pipe, eq. (41) reads 


.ald 


P; — Pp, = “F (note absence of g) 


Here the context shows that pressure is 
expressed in the dynamic unit |b. /ft.hr.” 
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All. Author All in discussing the fine- 


structure constant » 
h = Planck’s constant, e¢ = charge of electron) 
says: “it is frequently referred to as “‘ dimen- 
sionless,” for example when EppINGTON states 
that when we say that its value is 137 “ we are 
eliminating all reference to our artificial stan- 
dards.”’ It is clear, however, that this is not so.” 

We should like to say that the fine-structure 
constant, being the ratio of two like quantities 
(the velocity of light and the velocity of the 
electron in the first orbit in a hydrogen atom) is 
truly dimensionless so that its value is independent 


ch 
= velocity of light, 
~ 3 (c = velocity of ligh 


of our standards. 
It is only the expression for it that does depend 
ch 
on our standards (electrostatic egs system: on 
7 
,; 2ch 
practical system : “). 
* e 


The fact that the fine-structure constant is 
dimensionless should not be checked by sub- 
stituting in the e.s. cgs definitional equation the 
dimensional formulae for the various quantities 
from another system, involving ap ¢€), as it has 


been done. 


A 12. 
Pout [4]. 

Any charge can be expressed in e.s. cgs units 
or in e.m. egs units and it is therefore possible 
to equate the two, to give e.g. 


An interesting paradox is mentioned by 


1 e.m. egs unit of charge = 8 x 10"° e.s. egs 
units of charge 


(1) 


Now the e.m. unit and the e.s. unit are related 
to fundamental units as follows : 


(2) 


1 e.s. egs unit of charge = 1 cm*/* g'/? sec™ (3) 


1 e.m. egs unit of charge = 1 em'/? g'/? 


Inserting (2) and (3) in (1) gives: 


lsec. = 3 x 10'° cm (4) 


VERMEULEN [5] says the cause is to be sought 
in the essentially different symbolic nature of 
these products. 

Following statement elaborates hereon: Eq. 
(2) holds within the dimensional system used in 


forming e.m. units and within that system the 
product formations and so on can be regarded as 
normal mathematical routine. 

Eq. (8) similarly holds within the dimensional 
system used in forming e.s. units. 

In view of these restrictions the eq. (2) and (3) 
are in general incompatible and cannot be com- 
bined to give eq. (4). 

The e.s. and e.m. systems would become 
identical if such units of length and time were 
used that the velocity of light, c = 1 (see Part I, 
page 134). This is precisely what is found in 
eq. (4). 


A 18. Authors A 13 say that the diameter D, 
of bubbles formed in boiling is proportional to 


J 2 foo 
& (pi — Pv) 


and their list of symbols says : 
D, = diameter of bubble, ft. 
g = acceleration of gravity. 
£9 = conversion factor 
o = surface tension, dynes /cm or lb./ft. 
Pp» Py = densities of liquid and vapour, lb./cu.ft. 


If lb.,, and lb. are not distinguished, g and 
Z are identical. 


B. Inconsistencies in the substitution of numerical 
values in equations. 


B1. Authors B1 in studying the resistance 
? 
of pipes calculate a group 


rd and plot it 


2 
against a Reynolds number. 

The units for AP and p being derived from 
lb, and lb, respectively the above group should 
have been made dimensionless by including a 

Ib.,,, ft. 
factor g, = 82-2 __™ _ . 
a Ib., sec? 

The “friction factor’’- Re graph as given, 
accordingly differs from those obtained by a 
consequent application of coherent units. 

See also example B 6. 


B2. Authors B2 define the coefficient of 
friction y¥ for the turbulent flow of an incom- 
pressible fluid in a smooth pipe by 
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l 
Pi —P2=¥5. bpm 


== pressure drop (in cm water) over 
the length / (in em) 

d = diameter in cm 

w = mean velocity in cm/sec. 

p = density in g/cm*® 

ys = dimensionless. 


with p, — Pe 


With the other units as prescribed, p, — p. 
should have been expressed in dynes/cm* = g 
em™ sec.-*. The ratio between a pressure differ- 
ence and a liquid head is pg, i.e. for water, 981. 

B3. Authors B3 describe the addivity of 
resistances in mass transfer between phases. 
With the units given, k; is in cm/hr. and the 
product of the Henry coefficient with Ag is in 
litre /hr. em’. 

B 4. Authors B 4 give the following equation 
for the heat transport through a very dilute gas 
between parallel walls at temperatures 7, and 


T,: 
_8 cm _ 
@= 3p, [aR (7, -7,) 


where : p = pressure 
T = average temperature 
M = molecular weight 
R = gas constant. 

Q is not defined but the general list of symbols 
says: Q = quantity of heat, with dimension : 
keal. 

The Q in the above equation represents heat 
per area time but'Q and R can only be expressed 
in mechanical units. 


B5. Author B 5 gives for water at 100°C : 
specific heat 
¢, = 1-011 keal kg °C" 
viscosity 
p = 8-0 x 10° kg hr. m™ (i.e. kg, hr. m™*) 


thermal conductivity . 
k = 0-588 keal hr.' m= °C" 


and N,, (i.e. Prandtl group) 
= 1-75 (dimensionless). 
The Prandt! group had been previously defined 


as ““*. The above data would give ne = 


175 kgm 
1:80 x 10° kg, hr.2 
c, is thought to be referred to kg, or kg,,, the 
difference being that the forgotten multiplier is 
to be interpreted as g or g.. 

In order to arrive at the correct result : either 
the dimensional system should be adjusted to the 
units by the inclusion of the dimensional constant 
g or g., or the units adjusted to a dynamical 
dimensional system by expressing » in kg, m™ 
hr. 

Remark—This case is closely related to B11 
and B 12. 


B 6. 
by 





It is immaterial whether the 


Author B 6 defines the Reynolds group 


where : D = diameter of pipe (inches) 

u = velocity (ft. /sec.) 

$ = specific gravity referred to water at 
60°F (dimensionless) 

z = absolute viscosity in centipoises. 


The value differs from the usual one by a 


factor 


1 one 
~y The friction factor versus Re graph 


changes accordingly. 

It is felt that by this procedure the great 
advantage of dimensionless groups, their univer- 
sality, is lost. 


B7 and B8. Authors B7 and author B8 
prescribe to express diameter in ft., velocity in 
ft./sec., density in lb./cu.ft. and viscosity in 
lb./sec.hr. The Reynolds group thus calculated 
would be a factor 3600 low. 


B9. Authors B9 = B7 in their Notation 
prescribe noncoherent units, e.g. : 


particle diameter 
gas velocity 
density 

viscosity, 
diffusivity 


D,, mm 


pn,’ centipoise 
D,, sq.ft. br“ 
These cannot have served in calculating the 


dimensionless groups Re = D, Up and Sc = a 
p 
the numerical values of which are given correctly, 
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B10. Authors B10 in their Notation leave 
alternative choices, viz. diffusivity D in sq.cm/sec. 
or sq.cm/day; viscosity ms in centipoises or 


Ib. /sec.ft. The density p is expressed in lb. /cu.ft. 
It is also impossible to select a coherent set of 


data to substitute in £ 
pD 


B11. Author B 11 defines Pr = - (= 2) and 


= the Schmidt group. 


gives following data in his tables : 
= in keal/kg °C 
= in kg sec./m? 

in keal/m hr. °C 
= dimensionless 


; , . . 
However, with these units 7 is expressed in 


ec, 


hr. sec. pai , s 
we: and the multiplication with 3600 ; 


< 9-81 
P. 


m . 
——, is left to the reader. 
sec. 


B12. Author B12 does the same as author 


Bll. 


B18. Author B13 defines Pr = ¥ with the 
following units : 


c, in BTU/Ib. mole °F 
p in lb. /sec.ft. 
k in BTU/hr. ft. °F 


Note inconsistencies : lb. mole and lb., hr. and 
sec. 
The Se group cannot be calculated either. 


5. Tue SuITABILITY OF VARIOUS 
Systems or UNITS TO THE 
CHEMICAL ENGINEER 


In this section we shall: 


Enumerate the properties, desirable in a 
system of units ; 


Discuss these desiderata and examine in how 
far they can be reached ; 


Compare various systems. 


A. Desiderata 


For convenience and clarity the following is 
considered desirable : 


General requirements 

1. To use a coherent system of units in the parti- 
cular field in which the chemical engineer is 
working. 


2. To let such a coherent system of units cover at 
least the fields of mechanics and heat. 


. To strive to a one-to-one correspondence be- 
tween the quantities and their units and/or 
dimensional formulae, i.e. 

8a. Each quantity should have a single 
dimensional formula and unit. 

3b. Each dimensional formula and _ unit 
should correspond to a single quantity 
only (see FLEIscHMANN [7)). 

3c. As a logical consequence of 3b physical 
quantities should not be dimensionless. 


. To use a small number of dimensional constants. 


5. To employ such dimensional formulae and to 
express the derived units in such a manner that 
clarity is best served. This means that it should 
be obvious how the chosen expression is related 
to a definitional equation for the quantity 
under consideration. In this way it is easiest 
to remember the unit for the given quantity 
and vice versa. 

. To choose such a system that the acceleration 
by gravity appears only when gravity actually 
plays a role. 


Specific requirements re units 
7. Units should be well defined. 


8. Use must be made of following units : 
metric British 
length morem ft. 
mass and/or force kg or g Ib. 
time sec. or hr. 
temperature °C °F or °C 
. Due consideration must be given to the systems 
used by physicists, mechanical engineers and 
electrical engineers, since their equations and 


data will be used. 


B. Discussion of desiderata 
1. Coherent units 
A great majority of physicists and engineers use 
a coherent system of units. If their basic data 


sec. 
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are given in other units, they carry out the con- 
versions before substituting them in the equations. 
It should however, be placed on record that it 
is not necessary to do so (see discussion on THORNE 
and Watsnaw’s method, Part II, page 169). 


2. Width of field 

Special systems are sometimes used in a restricted 
field. Permeability of porous media, for instance, 
is expressed in darcy’s. We might say that there 


is an “atm./em; cm/sec.; centipoise; darcy 
system” for the use in Darcy’s Law only. 

As long as this is only applied in this particular 
field it works satisfactorily. If, however, per- 
meability is related to particle size, shape and 
porosity it becomes more practical to express it 


in m® or em®*. 


Likewise, in heat transfer there is a preference 
for hours, in other fields for seconds. Inconsis- 
tencies in the use of seconds and hours become 
apparent in examples B 7-8-9-11-12 and 13. It is 
considered that a consistent system for mechanics 
and heat is to be preferred and that it should be 
based on the second. 

The various static, dynamic and engineering 
systems, amplified with a unit of temperature and 
in some cases with a heat unit, as discussed, all 
fulfil this requirement. 

In the following it is taken for granted that 
such a coherent system will be used. 


3. “One-to-one” correspondence, and 


4. Small number of dimensional constants 


Rule 3a asks for a large number of dimensional 
constants ; a compromise with rule 4 must thus 
be found. 
Notable violations of 3a are : 
egs system applied to electricity and mag- 
netism, giving two sets of units ; 
English Engineering system (old and new), 
giving alternative choices for various units. 


Notable violations of 34 are that the same unit 


is employed 


wn for 
cgs e.s. system electric capacity and 


length 


in for 
egs e.m. system inductance and length 
English Engin. system, 


old version force and mass 


Notable violations of 3c are : 
egs €.s. system dielectric constant 
cgs e.m. system magnetic permeability 
heat units described in 

Part I, p. 138-9, ad 3 specific heat, specific 


entropy 


The necessity to use a small number of dimen- 
sional constants can be underlined by reference 
to frequent conversions requiring the mechanical 
equivalents of heat and electrical energy and to a 
considerable amount of misunderstanding about 
the conversion factor g, (Par. 4, Ex. A 1, 3, 4, 13). 

The practical system of course introduces jy 
and «, but avoids the conversion from practical 
electrical units to cgs units and back. 

Further refinement in the sense of requirement 
8 would be possible by distinguishing the three 
length co-ordinates in dimensional analysis and 
in units, or perhaps by introducing vector con- 
cepts. This is considered too complicated for 
general use, i.e. it is considered inevitable that 
torque and work have the same unit and that an 


angle (= arc/radius) will be dimensionless. 


5. Clarity of dimensional formulae 


Regarding clarity of dimensional formulae it is 
essential to introduce one new fundamental unit 
when entering a new field of physics, i.e.: 
temperature in the field of heat (all systems 
discussed, see Part I, p. 138) ; 
one electrical unit in the field of electricity 
and magnetism (practical system, see Part 
I, p. 139). 


This is contrary to the principle of making 
so-called based on 
alone, which for a long time has been the ultimate 


absolute units, mechanics 


aim of physicists. 


Requirements 3 and 5 are closely related, both 
acting contrarily to No. 4. 


6. The place of “g’’ in equations 


Requirement 6 disallows static systems (for the 
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reason explained in Part I, pp. 137-8) and the old 
version of the English Engineering system (Part I, 
p. 136). 

7. Definitions of units 

In order to properly define a unit of force, with 
the aid of a standard kilogram, or pound, it is 
necessary to specify a standard acceleration of 
gravity (g = 980-665 cm/sec.” at 45° latitude at 
sea level). Consequently the force of gravity 
acting on IIb... elsewhere in general is not 
1 lby. This reduces the supposed simplicity of 
static systems. 


8. Use of existing units 


All systems discussed in this paper conform to 
this requirement. 


9. Units used in physics and various branches of 
engineering 

In general physicists prefer dynamic systems and 

mechanical engineers static systems. 

Static systems have the advantage of having 
well understood units for force, 
(derived from kg,, gy, lb..) but the disadvantage 
of not having suitable mass units. This condemns 
them for use in general physics. 


pressure, etc. 


Dynamic systems have suitable mass units but 
the unit of force (newton, dyne, poundal) is con- 
sidered less suited by mechanical engineers. 

The English Engineering system allows con- 
ventional units for both force and mass, however, 
at the expense of conformity to requirement 8a. 


C. Comparison of systems 
It is considered that the English Engineering 
system (new version) has more advantages than 
disadvantages when compared to static systems 
and to the English Engineering system (old 
version). 
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Hence the choice is one between the English 
Engineering system (new version) — or its metric 
equivalent — and dynamic systems. 

It should be noted that E. Scumipr declares 
himself in favour of a kg,;, kg,,, m, sec. system 
but does not use it since there is no further 
literature in it. 

The comparison between engineering systems 
(mew version) and dynamic systems is different 
for British and metric units. 

In British units the argument to be in line 
with physics does not count so much, much of 
the physical data being in metric units. The 


entire field of electricity is in the metric system. 
Therefore the English Engineering system may be 
a very suitable compromise to satisfy the mech- 
anical engineers. 

Its disad vantages were noted to be the following : 


1. Alternative possibilities of defining dimensions 
of physical quantities with consequent alterna- 
tives where to place the dimensional! constant g.. 

2. Difficulty to grasp the significance of g.. 


A variety of false statements have been found 
(examples A 2, 3 and 4). 

It is necessary that the pounds (mass and force) 
are well distinguished. Simplest, but not self- 
evident, are lb. and Lb. [6]. In this paper, Ib, 
and lb, have been used. The abbreviations Ib 
(mass) and lb (force) are also employed. 

In metric units contact with physics counts 
much more. Here a dynamical system is consi- 
dered to be at its place. The choice between cgs 
and practical units in physics is largely governed 
by the fact that practical units are indispensable 
in electricity and signs are that cgs units will be 
completely abolished in this field. There is accord- 
ingly a strong current to do the same through all 
physics. Hence it is considered that the metric 
chemical engineer should use the practical system. 
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Mixing by agitation of miscible liquids 


Part I. 


J. G. van pe VusseE 


Koninklijke /Shell-Laboratorium, Amsterdam 


( Received 9 July 1954) 


Summary—tThe performance of several types of agitators has been studied for the mixing of 
miscible liquids in batch operation. 

The time was measured necessary to homogenize two super-imposed layers of miscible liquids. 
An optical technique, known as the “ Schlieren "' method, was used to determine the moment 
at which uniformity is reached. In addition, the power consumption of the stirrers was 
measured. The influence of liquid properties and of the conditions of agitation was ascertained 
and could be explained. 


The results enabled relations containing dimensionless groups to be derived for both mixing 
time and power consumption. 

An important result of this study is the discovery that in the region of turbulent flow the 
mixing time is related to the pumping capacity of the stirrer. The mixing time appears to be 
approximately proportional to the time required to circulate the liquid once, this time being 
equal to the ratio of the liquid volume to the pumping capacity of the stirrer. 


General rules for the optimum shape and dimensions of the stirrers were established from 
the relations derived for mixing time and power consumption. 


A few experiments were carried out on steady-state mixing (continuous mixing) of liquids, 
the results of which were in agreement with those of the batch experiments. 


Résumé—L auteur étudie l'efficacité de nombreux types d'agitateurs pour le mélange de liquides 
miscibles dans une opération discontinue. Il mesure le temps nécessaire &4 homogénéiser les 
deux couches superposées. I utilise une technique optique, connue comme méthode de 
“ Schlieren,” pour déterminer le moment of lhomogénéité est atteinte. I] mesure aussi la 
consommation en puissance des agitateurs. I] établit et peut expliquer influence des propriétés 
du liquide et des conditions d'agitation. 

Les résultats permettent d’établir des relations, impliquant des groupes sans dimension a la 
fois pour le temps de mélange et la puissance consommée. 

Cette étude apporte un résultat important : dans la région d’écoulement turbulent, le temps 
de mélange est relié & la capacité de pompage de l'agitateur. Le temps de mélange sembl- 
approximativement proportionnel au temps nécessaire & un seul cycle de circulation pour le 
liquide, ce temps étant étal au rapport du volume liquide a la capacité de pompage de l’agitateur. 

L'auteur établit des lois générales pour choisir la meilleure forme et les meilleures dimensions 
des agitateurs 4 partir des relation dérivées pour le temps de mélange et la puissance consommeée. 

Il a fait quelques expériences dans le cas, a I'état stationnaire, du mélange de deux liquides 
alimentant en continu une capacité, Les résultats obtenus concordent avec ceux des expériences 
en discontinu. 
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1. INTRODUCTION 
1.1. Definitions 


“* Mixing ” is understood to be any operation used 
to change a non-uniform system into a uniform 
one. A quantity of matter may be called uniform 
or homogeneous when the composition of a volume 
element of appropriate size does not deviate by 
more than a fixed amount from the average com- 
position of the entire system. 

Mixing is very often part of a chemical or 
physical process, such as blending, dissolvmg, 
emulsification, heat transfer, chemical reactions, 
etc. Chemical engineering recognizes it as one of 
the “ unit operations ”’ and it enjoys an increasing 
interest in literature. In spite of its importance 
mixing is still an “art” in most cases and the 
design of agitators and mixing vessels is largely 
based on experience [3] [6]. 

The devices for mixing are numerous ; they are 
greatly dependent on the purpose in view (e.g. 
heat transfer) and on the nature of the substances 


a PADOLE STIRRER bd. PROPELLER STIARER 


1. 
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to be mixed (state of phase, consistency). As to 
the latter point, the state of phase greatly in- 
fluences the choice of mixing method. Mixing can 
be classified according to the various combinations 
between the phases gas, liquid, and solid. Liquid- 
liquid mixing, for instance, is understood to be 
the mixing of two miscible or immiscible liquids. 
Important methods of mixing are: flow mixing, 
e.g. circulation by pumping, injection ; vibrational 
mixing, e.g. by ultrasonics ; mixing by rotating 
agitators [10]. 

In mixing two phenomena can be distinguished 
a splitting up of the volume elements and an 
interchange of these. In several cases, e.g. in 
mixing miscible liquids, the molecular diffusion 
also contributes to the ultimate homogeneity of 
the system to be mixed. In practice, the exchange 
by molecular diffusion, is negligible as compared 
with the eddy diffusion, due to the turbulence in 
the agitated system; it is only important when 
considering volume elements of molecular dimen- 
sions. 
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Mixing by agitation of miscible liquids—Part I 


1.2. Various types of agitators and their 
flow pattern 


Mixing by agitation of liquids involves the transfer 
of momentum from the moving stirrers to the 
liquid. According to the way in which this occurs, 
stirrers may be divided into two categories : 


(a) The momentum is transferred by shearing 
stresses, i.e. the transfer is perpendicular to 
the direction of flow. This category includes 
the cone stirrer, the bulb agitator, the rotating 
dise and others. 


The momentum is transferred by normal 
stresses, i.e. the transfer is in the direction of 
flow. This category includes the paddle 
stirrer, the turbomixer and the propeller 


(Fig. 1). 


The latter three may in their turn be regarded 
as representatives of types distinguished according 
to flow pattern. We may distinguish three prin- 
cipal possibilities of flow according to the three 
axes of the coordinate system : 

Tangential flow—This is found with all normal 
rotating stirrers, if no baffles or deflecting blades 
are present and the stirrer is placed centrally. 

Axial flow 
agitators. 


This is found with propeller type 


Radial flow—This is found with turbine type 
stirrers (often fitted with deflecting blades), and 
to a certain extent with paddle stirrers. 

As tangential flow does not contribute to the 
mixing in a vertical direction this type of flow is 
not particularly desirable. It will be shown later 
that tangential flow is always converted into 
radial flow by the centrifugal forces. This con- 
version can be promoted by baffling. In this way 
a “ top-to-bottom "’ mixing is attained, as in the 
case of the axial and the radial flow pattern. 


1.3. Scope of the investigation 


The investigations on mixing are always based on 


a mixing criterion which permits an evaluation 
As a rule, the choice of 


of the stirring process. 
criteria depends on the systems to be mixed 


together with the purpose of the operation. 

The present investigation is concerned with the 
mixing of miscible liquids. For this system the 
mixing time has proved to be an attractive 
criterion. The mixing time is the time required to 
produce a uniform mirture of two miscible liquids. In 
addition to the mixing time, the power consump- 
tion of the stirrer under various stirring conditions 
is measured, By comparing the mixing times and 
of the 


conclusions can be drawn as to their efficiencies. 


power consumptions various agitators 



































LIGHT SOURCE 


FIRST “SCHLIEREN” LENS 
SECOND “SCHLIEREN® LENS 


MIXING VESSEL WITH STIRRER 
SCREEN IN PLANE OF IMAGE OF V 


FIRST KNIFE - COGE 


INTERCEPTING HALF OF THE BEAM 


SECOND KNIFE - EDGE ( IN PLANE OF IMAGE OF M™,), 
INTERCEPTING THE OTHER HALF OF THE BCA! 


Fic. 2. Experimental assembly. 
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2. APPARATUS AND MEASURING 
TECHNIQUE 


2.1. Mixing time 


The mixing time was determined by what is known 
as the “ Schlieren ’’ method. The arrangement, 
which is a modification of that used by TorPLer 
[17], is represented schematically in Fig. 2. In 
order to avoid refraction the cylindrical mixing 
vessel was placed in a vessel with plane-parallel 
walls. Optical inhomogeneities in the liquid in 
the mixing vessel, in the form of gradients in 
refractive indices, will produce ‘“ Schlieren,” 
whereas absence of the latter indicates homo- 
geneity of the liquid mixture. 

In these experiments we started from a condi- 
tion in which the two liquids formed two super- 
imposed layers. If the liquids are carefully intro- 
duced into the vessel the interface is clearly 
visible. The following liquids were used : 
water and dilute solutions of acetic acid in 

water ; 


water and dilute solutions of glycerol in water. 


The mixing time was the time elapsing between 
the moment the stirrer starts mixing and the 
moment of disappearance of the “ Schlieren.” As 


to the reliability of the “ Schlieren” method for 
this determination, the following may be noted. 

It will be clear that the refractive index or the 
difference in refractive index has no influence on 
the time required to homogenize the two liquids. 
As stated above the “ Schlieren *’ method is based 
on the fact that rays of light are deflected owing 
to refractive index gradients in the mixing 
vessel. It is therefore conceivable that the 
difference in refractive index influences the 
determination of the mixing time. Experiments 
have shown, however, that only the “ intensity ” 
of the “ Schlieren ”’ is influenced by this difference. 

The mixing time as determined with this 
method proves to be independent of the difference 
in refractive index. Measurements made at a 
given difference in density give the same mixing 
time irrespective of the difference in refractive 
index. This phenomenon is in agreement with 
the observation that the intensity of the striae 
does not vary greatly during the mixing process 


until there is a sudden decrease at the end of 
the experiment (see figure). 


INTENSITY 
OF STRIAE 
(4M) 
(429),>(O9p)2 
( 


(2) 











|, MIXING TIME TIME OF M IxinG 


Fic. A 


This behaviour can be explained by a study of 
the mechanism of the mixing process. It will be 
shown later that the mixing time is more or less 
equal to the time required to circulate the liquid 
once under turbulent conditions. During that 
time a continuous flow of liquid passes the stirrer, 
giving a striae image of fairly constant intensity. 
When all the liquid has passed the intensity 
decreases abruptly. 

The moment of disappearance of the striae is 
dependent on the sensitivity of the method. How- 
ever, if the moment immediately after the sudden 
decrease in intensity of striae (shown in the figure) 
is taken as the end point, the influence of the 
sensitivity is very small. It appears that the 
reproducibility is lower than 5% for various 
sensitivities of the set-up. The apparatus detected 
striae formation occurring at differences in refrac- 
tive index of 10° under normal operating con- 
ditions. The difference in refractive index between 
the two starting liquids is of the order of 10°. 

For laminar flow, the end point of mixing is 
very difficult to determine, as the transition to the 
homogeneous state is a gradual one. This is closely 
connected with the mechanism of the mixing 
process. Two effects play a réle in reaching the 
uniform state: the long range convective mass 
transfer effected by the pumping action of the 
stirrer and the (eddy) diffusion between the small 
elements of volume. 

In the case of turbulent flow the diffusion is so 
great that the mixing time is governed by the 
pumping action of the stirrer. With laminar flow 
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the diffusion is so small that the mixing time is 
controlled by the diffusion effect. 


2.2. 


In the present investigation the following 
methods have been applied for the experimental 
determination of the power consumption : 


Power consumption 


(a) The stirrer is driven by a falling weight con- 
nected to the stirrer spindle by a piece of 
string. After release the weight first falls with 
an accelerated motion, subjecting the stirrer 
to a frictional torque in the liquid which 
increases with the speed and hence with the 
rate of fall. As soon as this frictional torque 
has become equal to the torque caused by the 
weight the motion will be uniform. The dis- 
tance travelled by the weight per second 
multiplied by the weight gives the power 
consumption of the stirrer. Naturally, the 
experiment must first be carried out with 
the stirrer rotating in air to determine the 
frictional loss in the bearings, so that a 
correction can be made for this loss. 


The mixing vessel is placed on a rotary table. 
If the axes of stirrer and rotary table are 
parallel the couple required to prevent 
rotation of the vessel is equal to the couple 
exerted by the stirrer (the whole system is 
at rest and the resultant is therefore equal to 
zero). The angular velocity multiplied by the 
couple is equal to the power consumption of 
the stirrer. This determination must also be 
corrected for friction in the bearings and is 
only applied in the case of relatively great 
torques. For small torques the accuracy is 
not sufficient. 


8. MIXING PERFORMANCE OF PADDLE 
STIRRERS, PROPELLER STIRRERS AND 
TURBOMIXERS 


3.1. Mixing time 


8.4.1. 


The mixing time (7') was determined as a function 
of a number of variables. It was found to depend 
on the following liquid properties : 


Dimensional considerations 


# = final viscosity of the mixed product 
p = density of mixed product 


Ap = difference in density between the liquids to 
be mixed and on the stirring conditions : 

n = number of revolutions of stirrer per unit 

time 
= diameter of stirrer 
D = diameter of vessel 
H = height of liquid in vessel. 

Moreover, it must depend on g = acceleration 

of gravity. 

Other factors which influence the mixing are 
location and mounting of the stirrer (eccentricity), 
the ratio by volume of the liquids to be mixed, 
the shape of the vessel (presence of baffles) and 
the shape of the stirrer. 

Let us first consider the mixing time T as a 
function of (n, uw, p, Ap, d, D, H, g). Together 
with 7 these quantities can be combined to 
9-3 = 6 dimensionless groups. For these groups 
we have taken : 


an where V = volume of the liquid 


(1) 
(V ~ D* H) 


(2) pnd = Reynolds’ number 


ms 
n*d? 


—. == Froude number 


(3) 
(4) 
(5) 


(6) 


The above quantities have been adapted to the 
mechanism of mixing, so that they also have a 
physical significance. It has already been said 
that the mixing time is mainly determined by the 
pumping capacity Q of the stirrer. Accordingly, 
the mixing time is proportional to the time 
required to pump all the liquid round once : 
7 
T~ 3 Thus, r = cb represents the number of 
times the volume of liquid has to be pumped 
round. Dimensionally, the expression for the 
pumping capacity has the form Q = nd*. How- 
ever, it is obvious that in correlations use must 
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be made of expressions for the theoretical pumping 
capacity §th which take into account, not only the 
diameter d of the stirrer, but also type, shape and 
other dimensions (see table, Pumping Capacity 
3.2.4). The dimensionless mixing time 7 is then 
given by —— 

In these expressions for the pumping capacity 
only stirrer dimensions are involved. During 
mixing, however, other factors also affect the 


actual pumping capacity. For example, the 
rotating stirrer has to overcome a certain counter- 
pressure, because liquid of lower density is pushed 
into liquid of higher density, thus reducing the 
pumping capacity. This influence is characterized 
by a density Froude number giving the ratio of 
the dynamic head of the flowing liquid, which is 
assumed to be proportional to pn*d*, to the static 
head ApgH that has to be overcome. This 
quantity corresponds to the product of the 


VARIATION OF &P AND ys axiom? 
VARIATION OF 4 AND w 


Va 14x10 me? 
v2 45x10 2m? 
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Fic. 8. Paddle stirrer with a right blade angle. Influence of the modified Froude number 
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dimensionless groups mentioned in (8) and (4). 
al 


With decreasing values for the quantity A pg 


the mixing time will naturally increase. 
The dimensionless mixing time 7 will also 
depend on the condition of flow in the mixing 


a ri 
= | » = ) 


vessel, i.e. on Reynolds’ number 


NUMBER OF VANES 
@ = 100 mm 
we 4 mm 
B= «6° 


Finally, various length ratios play a rdle, 
notably those which give the dimensions of the 
stirrer in comparison with the total liquid volume. 
If the vessel is too large or of unsuitable dimen- 
sions in relation to the dimensions of the stirrer, 
“dead corners’ with little turbulence will be 
formed, which lengthen the mixing time. 

The expression for the mixing time may now 


$:12 


f = 2 (FACTOR RELATED TO NUMBER OF VANES) 


x V=7x0 3m? 
O v= 14002 m? 
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Fic. 4. Turbomixer. Influence of the modified Froude number on the dimensionless mixing time. 
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nd pn*d d H nd* pn*d* d 
be formulated as =f ao : t=; (“") (ea) 4 (pares) (1) 
During the correlation of the data we tried to 
represent the right-hand side of the equation by It is obvious that this will not always be 


a product of functions, viz. possible, although here good results were obtained. 
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Fic. 5. Paddle stirrer with inclined blades. 
Influence of the modified Froude number on the dimensionless mixing time. 
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8.1.2. Influence of the density Froude number _ based were carried out in a region where Reynolds’ 
number and the length ratios have little influence. 
In a large region the relation may be represented 
This relation is given for the various types of by ® power function : 
stirrer in Figs. 3, 4, 5 and 6. 
The measurements on which the graph are 
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Fic. 6. Propeller'stirrer. Influence of the modified Froude number on the dimensionless mixing time. 
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where for a paddle stirrer 0-3 


a paddle stirrer with 
inclined blades 


a turbomixer 


a 
a 
a 


a propeller stirrer = — 0-25 


For Ap = 0 the mixing time 7 approaches a 
finite positive value, which may be determined 
by extrapolation to Ap = 0 in a T vs. Ap graph ; 


hence, for very high = 


fa values + approaches 
A pg ‘ 


H 
a limiting value. 

We have also applied an injection method in 
some cases. Mixing-time measurements were 
carried out by suddenly adding a small amount 
of a glycerol-water mixture to the water during 
stirring, and measuring the time necessary for the 
striae to disappear. The results are indicated in 
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Fic. 7. Paddle stirrer with right and with inclined blades. Influence of Reynolds number on mixing. 
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Figs. 3, 4, 5 and 6. The values found are invari- In addition, the influence of the volume ratio 
ably higher than the limiting values found by between the two liquids to be mixed was ascer- 
extrapolation to Ap = 0 in a T vs. Ap graph, tained. As will be demonstrated in section 3.6, 
which suggests that in our experiments the the increase in potential energy for a volume 


*t* 
quantity pn still has some influence. ratio of 7 £ J between the two liquids to be 
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Influence of Reynolds number on mixing. 
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mixed, at a constant total liquid volume, is equal 
to 4¢(1 — ¢) ApgH V. Now if for the density 


Froude number we tak = 

roude number we ta ‘ia — $) pg 
pn*d? es 
of A pel (which values are equal for the normal 
case of ¢ = 0-5), the results of the experiments 
with a varying ¢ can be correlated in the normal 
graph of dimensionless mixing time versus density 
Froude number. 


instead 


3.1.8. Influence of Reynolds’ number pnd 
pe 


This relation is given in Figs. 7 and 8. For values 


PADOLE STIRRER 


© VARIATION OF 4 (0 = 148mm H= 116 mm) 


& VARIATION OF 0 ANOH (d= 45mm) 


GO BADGER AND al, (CH ENG. 27 1176 ‘22)(1) 


ne 104-105 


Te? wee? wong?) 23 
Vv ss 


2 
Re > 10° to 10* the quantity +. (ee, 


to be nearly constant. 

For lower values it increases rapidly (7. Fr®® 
~ Re*). The type of flow is laminar, with per- 
sistent liquid eddies that show little exchange 
with their surroundings. Inhomogeneities are 
particularly persistent in the “ dead corners.” 

The behaviour of the streamlines also revealed, 
for all stirrers studied, a transition from the 
laminar to the turbulent state. The streamlines 
could be made visible by introducing small gas 
bubbles into the rotating liquid. At a certain 
critical stirrer speed the track of the gas bubbles 


PROPELLER STIRRER 


X 3.0L dup = 0x8.104 me 
VARIATION 0 AND H 
Re 104- 105 


prta? } 0.25 


2 
}) { apen 














ys 







































































TTT TTT 






































aati 
























































= : , 


ots 


Fic. 9. Paddle stirrer and propeller stirrer. Influence of the vessel dimensions on mixing. 
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becomes irregular. The change gradually set in 
at a value of the Reynolds’ number : 
=> 250. 

In the laminar region the mixing mechanism is 
no longer controlled by the pumping capacity but 
by the exchange between the more or less per- 
sistent liquid eddies. The concentration at a 
given point in the vessel approaches the equili- 
brium concentration much more gradually and 
there is no question of a sharp transition from 
not-yet-uniform to uniform composition of the 
liquid, as is the case in the turbulent range. In 
the laminar region sharp measurement of the 
mixing time from the striae image is not possible. 


Re scndtien 


3.1.4. Influence of dimensions of vessel relative to 
stirrer dimensions 


The relation is shown in Fig. 9. It was possible 


to include the and ratios in one combina- 


3 
tion, viz. p?* H°* This quantity may also be 


a 


y and more or less represents 


the ratio of stirrer volume to vessel volume. In 
3 


d 
= : ' 
the range 10° < DP 7s < 4 
2 
quantity +r. ( al 
A pgH 
8 
outside this range it increases. For DP* A < 10% 


H 05 
written as . 
D 


107 the 


~<a 
is practically constant ; 


the vessel volume is so large that dead corners are 


formed which lengthen the mixing time, whereas 
8 


d 
for ps yes > * x 10 the vessel volume is so 


small in relation to the stirrer volume that the 


movement of the liquid by the stirrer is impeded 
8 


by the vessel wall being so near. For p?* H®* 


> 10° the relation can be represented by a power 
function : 


pnd? -a 1 da +1 
*(Sin) ~ (ora) 
Influence of location of stirrer, etc. 


8.1.5. 


Eccentricity of stirrer—From the curves shown in 
Figs. 3-9 it is apparent that the behaviour of all 


types of stirrers is about the same. However, 
differences in the behaviour of stirrers occur when 
they are positioned eccentrically. 

By eccentricity we understand the distance 
between the stirrer shaft and the axis of the 
mixing vessel divided by half the diameter of the 
vessel. The influence of the eccentricity on the 
mixing time is shown in Fig. 10. Unlike the 
influences discussed above, this influence proves 
to be greatly dependent on the type of stirrer. In 
the case of stirrers causing axial flow, like the 
propeller stirrer and the paddle stirrer with 
inclined blades, eccentricity shortens the mixing 
time. This is due to an increase in the relative 
velocity of the rotating stirrer and the liquid 
rotating with it. With an eccentric stirrer this 
relative velocity will be greater than when it is 
centrally placed. As will be shown in section 
3.2.3 greater relative velocity causes a greater 
axial pumping capacity and therefore a shorter 
mixing time. In the case of stirrers causing tan- 
gential and radial flow, like the paddle stirrer 
(with straight blades) and the turbomixer, 
eccentricity lengthens mixing time. The tangen- 
tial and radial pumping capacities also increase 
as the relative velocity increases. As, however, 
here the vicinity of the wall impedes the flow of 
liquid in a radial direction the mixing time 
becomes longer. 

In the case of a paddle stirrer with the blades 
set at an angle « of 40° to 50°, the mixing time is 
found to be independent of the eccentricity, hence 
it may be assumed that here the two effects 
counterbalance each other. 

Influence of the clearance between stirrer and bottom 
of the vessel—Kramers and Know. [8] have 
ascertained this influence for a propeller stirrer. 
Their measurements show that a clearance equal 
to half the liquid height generally ensures the 
shortest mixing time. Our own experiments with 
a paddle stirrer pointed in the same direction 
(Fig. 11). The investigation into the influence of 
this variable was therefore discontinued. Accord- 
ingly, in all experiments the clearance was equal 
to half the liquid height, unless stated otherwise. 

Influence of inclination of the stirrer shaft—Rota- 
tion of the liquid with the stirrer can be sup- 
pressed, not only by an eccentric position of the 
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stirrer, but also by placing the stirrer in the 
mixing vessel in an inclined position. When the 
stirrer shaft of a propeller is set at an angle of 
30° or more to the axis of the mixing vessel the 
mixing time becomes about 0-7 times that under 
normal conditions. This effect, therefore is of 
the same magnitude as the influence of eccentricity 
of the stirrer. 

Influence of baffles—In the experiments with 
baffles only the assembly shown in Fig. 14 has 
been studied. This type of baffles (mounted at 
right angles to the wall of the vessel) is commonly 
used in practice. 


The object of placing baffles in the mixing 
vessel is to prevent the liquid from rotating with 
the stirrer. The ensuing advantage is that this 
counteracts vortex formation, which beats air 
into the liquid. Moreover, baffles cause the flow 
of liquid to be directed more axially, which 
generally has a favourable effect on the rate of 
mixing. In the case of propeller stirrers baffling 
has the same effect on the mixing time as eccen- 
tricity of the stirrer, especially at higher rotor 
speeds. The relative velocity of the stirrer and 
the liquid then approaches the stirrer speed, 
resulting in a higher pumping capacity. At lower 
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rotor speeds baffling has less influence on the 
flow pattern and on the mixing time. 

With paddle stirrers and turbomixers baffling 
has about the same influence on the mixing time 
as with propeller stirrers. Turbomixers are often 
provided with deflecting stator blades, resulting 
in a perfectly radial flow pattern. 


8.2. 
(influence of shape and dimensions of stirrer) 


3.2.1. Paddle stirrers and turbomizers 
When calculating the dimensionless mixing time 


Pumping capacity 


TQth 

a 

must be expressed in characteristic and measur- 
able quantities of the stirrer. 

We will successively deal with the pumping 
capacities of the paddle stirrers (with straight 
blades), the turbomixer, the paddle stirrer with 
inclined blades and the propeller stirrer. 

The paddle stirrer with blades set at right angles 
imparts to the liquid a momentum that gives it a 
tangentially rotating movement. From the point 
of view of mixing this movement is undesirable, 


the theoretical pumping capacity, Qth, 


R= CLEARANCE ( DISTANCE STIARER TO BOTTOM) 
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Fig. 11. 


Influence of clearance h on mixing time. 
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because it does not contribute to mixing in a 
vertical direction, and hence cannot counteract 
the influence of gravity. The fact that the 
paddle stirrer nevertheless has a good mixing 
effect is due to a secondary radial flow pattern 
formed on account of the centrifugal forces set 
up. 

The flow pattern of the turbomixer greatly 
resembles that of the paddle stirrer. Here, too, 
a radial flow is set up by centrifugal forces. Only 
the pumping capacity due to this radial flow has 
been considered. Since shape and action of the 
turbomixer are identical to those of the centrifugal 
pump the expression for the pumping capacity of 
a centrifugal pump can be used. 

From the theory of centrifugal pumps it follows 
that the so-called “ Euler’s head ”’ of the liquid 
leaving the centrifugal pump is 

__ Cg Wy COS Hy 


h ie 
_ g 


where c, = absolute discharge velocity of the 
liquid from the vane, 
a = the angle that c, forms with the tan- 
gent to the circle having a radius r,, 
r, = radius of vane wheel (see figure). 


equal to 


. 2p 


w = ANGULAR VELOCITY AS Eh 


@ = VANE ANGLE 


Fic. B. Discharge velocity diagram of turbomixer. 


According to BerNouuui’s law: hg = hy, + 
hat if friction losses in the pump are neglected. 

When no deflecting blades are present 
hyin = 4¢,"/g. The static head h,,, is assumed 
to be equal to the head }(w,/r,)*/g, if w, = 
the angular velocity of the rotating liquid in 
the vessel. 


Hence, pty wr COS ag = § pty” + $ p (wzr3)°- 


If == = q(q < 1), it follows that v,., = c, sin 
t = wr, sin 8 4/1 — g®, where 8 stands for the 
vane angle between the blades and the tangent 
to the circle having a radius r,. When q = 0, 
% = y and c,sin a, = wr, sin B. 

The radial pumping capacity is the product of 
the annular area 27r,w and v,,,4 


where w = width of stirrer blade ; 


therefore the radial pumping 
4n* nr? w sin B 4/1 —¢ 7 nd* w sin B 
/1—q*. For a paddle stirrer B = 90° 
sin 8 = 1, therefore Q = 7° nd* w 4/1 — q?. 
If all the liquid in the veses] is rotating at an 
1 and Q = 0. 
0 and Q = n* d* wn 
TQth 
V 


capacity = 


and 


angular velocity w, q = 
If baffles are present q¢ = 


sin 8. The dimensionless mixing time + = 


Tnd* wn’ sin B 
y 

Tnd* wn* 

ae 
When deflecting blades are present (as is often 
the case with turbomixers), again hg - 
how If v, is the discharge velocity of the liquid 
from the stator, h,;, = 4 v,". If no rotation of 
the liquid in the occurs hy, may be 
assumed to be zero, hence pc, wr, Cos a = }pv,.. 
—ke,sina, if k=—*-. The radial 

stator 

becomes 


becomes 7 = and in the case of 


paddle stirrers r+ = 
Nyin 5. is 


vessel 


Now UL 
velocity c,sina now Cc, sin a = 
wt, (—1 + 4/1 + 2h? tg? B) 


k? tg B 
0-7, so k? = 


In many construc- 


tions k = 4, which gives: 


cg sina & 2 ar (: 


Compared with the turbomixers without de- 
flecting blades, the pumping capacity is now 


larger by a factor Theoretically, the 


2 
1 + cos B 
pumping capacity would become infinitely large 
for 8 = 180°. In practice the velocities for 
B > 90° become so high that friction losses become 
important. As $ pc,” + (Agtat) = } pv,’ + 0 
and ¢? > vy) (har)... <9, which means that 
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liquid will flow into the slit between rotor and an are) revealed that including a factor sin 8 in 
stator. Consequently, (hy,,) will increase, the formula for the dimensionless mixing time 


rotor . . ““ a " 
and v,, and hence the pumping capacity will also has experimental significance (Fig. 12). The 
l 
, Suggests 
1 ve 

linearity between sin 8 and Q. Extrapolation to 
sin 8 = 0 produces a finite value for 7, although 


decrease. linear relation between sin 8 and 
A gradual transition from paddle stirrer to 
turbomixer (without deflecting blades) could be 
studied by curving the blades and investigating 
the influence of 8. 
Experiments in which the influence of the vane 


the theoretical radial pumping capacity would 
then have to become zero. Apparently, there is 


angle 8 of a two-bladed stirrer on the mixing time still some flow which causes mixing. 


was determined (the vane blade being curved in The influence of the blade width w on +r was 


$ Ww sec ~! 
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Fic. 12, Turbomixer 2 vanes. Influence of the curvature of the vane (vane angle 8) on the mixing time. 
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determined experimentally. It was found that Q 

, , u 

is proportional to w for values of ] 

( 

which is in agreement with the above-mentioned 

T@wnn*® .. w 
- For 


Vv d 


less than 0-5, 


expression 7 = values of more 


= es = ie , 
than 0-5, + increases with increasing i” Le. @ ts 
tf 


then less than proportional to w. 
3.2.2. Paddle stirrer with inclined blades 

A type of stirrer intermediate between the 
paddle stirrer and the propeller is the paddle 
stirrer with blades set at angles of less than 90°, 
This might also be called a propeller stirrer with 
a pitch decreasing towards the shaft. The relation 
between blade angle « and pitch p is given by : 


where r is the distance from the stirrer shaft. 

This type of stirrer produces both an axial and 
a radial flow pattern. Depending on the blade 
angle one flow pattern will predominate. Gen- 
erally, both occur and there is thus an axial and 
Q,, and Q.,. The 
total pumping capacity Q = Q,, + Qi The 
1 l l 

T Tox Trad 

if T,, is the mixing time when only axial flow 
occurs, and T',, the mixing time when only radial 
flow occurs. 

The influence of the blade angle « on the axial 
and radial (tangential) pumping capacity can be 
determined (theoretically) from the figure below. 


a radial pumping capacity : 


mixing time 7' can be found from 


a= BLADE ANGLE 
ws BLADE WIDTH 
v= VELOCITY 


Ta cos? a 


Velocity diagram for a paddle stirrer with inclined blades. 


The relative velocity v, between paddle and 
liquid is equal to the paddle velocity v, minus 
velocity v, of the liquid rotating with the stirrer 
(assuming that v, has no axial or radial com- 
ponent). Actually the liquid velocity v, will form 
an angle with the paddle velocity v,. The relative 
velocity wv, then becomes the vector difference 
v, — v, and can be split up again into the com- 
ponents perpendicular to and parallel to the 
paddle blade. It is now assumed that the momen- 
tum of the component perpendicular to the paddle 
blade is destroyed. The result is then a velocity 
of the liquid relative to the paddle in the direction 
of the blade, which is equal to v, cos a. The axial 
component becomes v, cos «sina and the tan- 
gential component v, cos? x. The absolute liquid 
velocity in tangential direction is then 


" , ens* 7 cin? | 
Up — U, COS” & vU,sin' a + 0, 


The tangential velocity is reconverted into a 
radial velocity by centrifugal action. 

Since the total liquid mass already has a 
the tangential speed 
relative to the surrounding liquid becomes equal 
to v, sin? « and the radial velocity, which is equal 
to the relative tangential speed, also becomes 
The relative velocity v, still 
depends on various factors. With an eccentric 
stirrer position and in the presence of baffles, v, 


tangential velocity  v,, 


equal to v, sin® a. 


will increase and become approximately equal to 
the paddle velocity, resulting in an increased 
pumping capacity. 

When the entire volume of liquid between the 
paddle blades has an axial velocity equal to v, sin 
a cos « and a radial velocity v,sin*« the axial 
pumping capacity becomes : 


id 
Q.. -| w, ($d) sin « cos « md dD ($d) 


0 


, 1 
= w, SiN & COS & == md® 


The radial pumping capacity becomes : 
Qua = = o, ($d) sin? « wdw sin « 
= $a, sin® « d*w 
For « = 90°, Q,, = 0 and Q,,4 = 2°n, dw. 
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Actually, it is not the entire volume of liquid 
” dw sin « that will have the above-mentioned 


velocity, but only a part y of that volume, which 
may be represented by the sector with an angle 


mw, sin « cos a d® 
» 


1 
of y x 360°. Then Q,, = y. ; 


and Q..4 = y~- $7, sin® « dw. The factor y 
will be influenced by the number of blades, and 
when this is small (1 to 3) it will be proportional 
to it. Now if we assume that the mixing time ts 


inversely proportional to the total pumping capa- 


© sP= 40 ag/a* 


x = I®JECTION METHOD 


4s 5-10%m 
w 2 075-0 em 
as 430 (pm 


city (Q,, + @,,4) and that the relative velocity 
v, is proportional to the circumferential velocity 
of the paddle, T becomes 


” . . 
nny d* sin «cos a + mn yd*w sin*® 2| 


Tt ~ (; 
o 
Experiments were carried out with a three- 
blade paddle, d being = 50 mm and w = 12 mm. 
The variation of T with « is then given by 
T ~ (0-7 sina cosa + sin*«)'. Experimental 
testing of this formula shows that the relation is 
valid for 0° < « < 60° (Fig. 13). In the range 
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Fic. 13. 3-Blade paddle stirrer with inclined blades. 
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influence of the blade angle « on the mixing time T. 
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60° < «a < 90° the experimental mixing time is 
shorter than that found with the aid of the above 
relation. 

Correlation of the mixing time with sin a shows 
that for 40° < « < 90° it is true that T ~ sin* « 
(Fig. 14). For smaller angles T ~ sin" «. 

The dimensionless mixing time 7 used in the 
correlations was finally defined as follows for 
paddle stirrers with inclined blades : 


Tnd*w n? sin *x 


T= p (k being a function of «) 


assigning to sin*x the value belonging to the 


angle « concerned (1 <k < 2). 
In the case of paddle stirrers with a blade angle 


90° the dimensionless mixing time is again 


Tr@wr* 


reduced to the value already found, 7 = —— y 


¢= 


3.2.3. 


The common propeller stirrer greatly resembles 
an ordinary ship’s screw. It has screw-shaped 
blades with a constant pitch p. If the propeller 
is regarded as a paddle stirrer with blades set at 
an angle we may here again derive an expression 
for the pumping capacity in accordance with what 
was found in 3.2.2. But here the blade angle « 
it varies in such a manner that 


Propeller stirrers 


is not constant : 


tan « = , r varying from zero (near the shaft) 


2nr 
to $d (extremity of stirrer blade). Now : 


3- BLADE PADOLE STIRRER 
dxuws Sx 0.75x0% at 


aPs= 40 kg/m? 
-3 3 


Vz 2.0 ° am 
as 430 rpm 







































































Fic. 14. Paddle stirrer with inclined blades. Influence of blade angle o on mixing time T. 
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P 


sin a = - 
Vv p® + 4n*)? 


2nr 
V/ p® + 4n%?2 
The expression for the axial pumping capacity 
would become : 
hd 
Qa. 


| 2rw sin acosar® dr 


5 


. In I 


L w ad’ f 
7 


For the radial pumping capacity the expression 
for the radial velocity is found from : 
r 


4 pv" | pw* sin® x r)r 


rad 


* 
0 


j pw" = In l 
tr- 


radial then becomes 


The 
Read 


pumping capacity 


Uraq TAw Sin & 
w dpu 2 
: - Jin( : 3) - 
7° a* 
2 /i ' a? 


Hence, in spite of the rather simple flow picture 
on which the calculation is based the formula for 


the pumping capacity has already become fairly 
complicated. Therefore, these formulae are not 
suitable expressions in the dimensionless mixing 
time of a propeller stirrer. Considering that the 
flow picture of a propeller stirrer is mainly axial 
and that it follows from the equation for Q,, that 
@. % Tnpd* 
ax V 
as the dimensionless mixing time. The pumping 
capacity of propeller stirrers was experimentally 
determined by Rusuron and co-workers [14]. They 
found that under optimum conditions the axial 
pumping capacity is directly proportional to the 


npd*, the expression r was chosen 


stirrer speed and to the square of the stirrer 
diameter. In later work they also found a 
proportionality between the axial pumping capa- 
city and the pitch of the propeller [15]. 

These results give also experimental strength 
to the assumption that Q = npd*. 


3.2.4. Recapitulation 
The expressions for the dimensionless mixing time 
in their relation to the stirrer type may be 
summarized as below. 

The first three formulae (paddle stirrers and 
turbomixers) may be classified under one general 


formula : 
Tnd® wr sin B sin *x 
\ 
which gives the mixing time for paddle stirrers 


(straight or inclined blades) as well as for turbo- 
mixers. This value of +r substituted in (2) gives a 





Type of stirrer 





Paddle (a = 90°) tangential 


Turbomixer (without deflecting blades) radial 


tangential 


Paddle (x < 90°) end extol 


Propeller axial 


Flow pattern 


Dimensionless mixing 
Nomencicture 


stirrer diameter 
width of blade 
vane angle 
— 
Tnd? w = sin B blade ungle 
} pitch 
Tnd? w n* sin *x 
—y 
Tnd* p 
y 
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general formula for the mixing time, which is 
valid in the turbulent region for a vessel of 
appropriate size : 
Tnd* wn’ sin B sin *x (p n*d* 
V (—— 
, T Ph 2 df? 0-25 
propeller stirrers in ; P (= » 9. 
I A pgh 
Number of blades of the stirrer—Generally, paddle 
and propeller stirrers have either two or three 





0-38 
) =~ 7, whereas for 


REYNOLOS 
at a2 


blades. In the case of turbomixers the number of 
blades (vanes) is usually much larger (12-20). The 
influence of the number of blades on the pumping 
capacity was not extensively investigated. It was 
observed that experimental results correlate well 
when the pumping capacity is taken proportional 
to the number of blades, at least, when the latter 
is small (1, 2 or 3). When it is large, as for the 
turbomixer, this proportionality no longer holds 
true; the pumping capacity then approaches a 
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Fic. 15. Vortex formation. Relation between Reynolds and Froude number at vortexing conditions. 
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maximum. The measurements with the turbo- 
mixer (12 vanes) may be made to coincide with 
the results obtained with the (two-blade) paddle 
stirrers if the pumping capacity is multipled by 
a factor 2 (f = 2). 


3.3. Vortex formation 


In most cases vortexing is undesirable because 
then gas is dispersed into the liquid. Moreover, 
mixing is not further improved when the speed 
of the stirrer n is increased beyond the critical 
speed n* at which vortex formation starts. Results 
of measurements on vortexing caused by centrally 
placed stirrers could be fitted into a correlation 
between a dimensionless group of Froude, viz. 
nd d_ n**@? pn*d* 
“g oh gh 

(Fig. 15) (h = distance between stirrer and level 


of liquid). For a Reynolds’ number > 10* the 
on 
uantity — 

q * A 


o2 
Reynolds’ values (Re > 10°) - = ~|( 
&g 


, and of Reynolds’ number 


For low 
pn*d* ) ° 
‘ 


proves to be constant. 


With the aid of these correlations the conditions 
to be used when vortexing must be avoided can 
be calculated in a very simple way. 


3.4. Rules for scaling up 


If a stirring apparatus is enlarged geometrically 
the conditions that have to be fulfilled to main 
tain the same mixing time in the large and smal! 
apparatus can be derived from (2). 
oa 10 = Leat\™ 
For Re > 10°-10*, pe ~ (Sn) (the ex 
ponent 0-3 being an average value) or 7 ~ n 
d°*, which means that (nd°*) must remain equal 
for T = constant. In practice this means that 
the number of revolutions of the stirrer should 
be nearly the same in the large and small appara 
tus. Generally, this claim will cause difficulties 
as the power consumption will then become too 
; ' one ind pn*d* y"" pnd*\~* 
high. For Re < 250, > ~(2 5 ( : 
or T ~n** d**, which means that (nd'*) must 
remain constant for 7 = constant. 
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